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ABSTRACT 


This  research  examines  the  applicability  of  a  simple 
objective  analysis  technique  on  a  personal  computer  to 
determine  the  detailed  distribution  of  surface  winds  during 
strong  downslope  windstorms.  A  three-dimensional,  mass-con¬ 
sistent,  diagnostic  wind  flow  model  applicaole  to  complex 
terrain  (COMPLEX)  was  modified  to  analyze  surface  wind 
fields  for  two  severe  chinook  windstorms  in  southern 
Alberta.  The  principal  modification  to  the  model  consisted 
of  specifying  the  vertical  structure  of  the  mesoscale  air 
flow  with  a  quasi-periodic  function  based  on  generally  known 
flow  characteristics  in  mountain  waves.  Specifically,  the 
boundary  layer  was  assumed  to  be  perturbed  by  a  two-dimen¬ 
sional  'hydrostatic  wave'.  Sensitivity  tests  were  performed 
to  determine  the  wave  length  and  boundary  layer  depth  which 
would  provide  the  best  estimates  of  observed  surface  wind 
speeds  at  selected  locations.  Results  of  the  most  success¬ 
ful  sensitivity  test  were  then  applied  to  a  sequence  of 
hourly  data  sets  for  each  windstorm  case.  Space-time  analj  - 
ses  of  observed  hourly  averaged  wind  speeds  were  implemented 
to  approximate  the  movement  of  the  hydrostatic  wave  in  each 
case . 

Realistic  results  were  produced  in  the  case  study  anal¬ 
yses  where  maximum  winds  followed  the  position  of  the 
'hydrostatic  wave'  trough,  with  the  highest  model  wind 
speeds  produced  when  the  wave  trough  was  closest  to  the 
mountains.  Problems  caused  by  the  assumed  two-dimensional 
wave  structure  included  unrealistic  changes  in  wind  direc¬ 
tion  at  some  boundary  points  where  the  model  flow  surfaces 
passed  beneath  the  terrain  surface.  Similarly,  the  inter¬ 
section  of  flow  surfaces  with  the  terrain  occasionally 
caused  wind  speed  maxima  to  occur  between,  rather  than  over, 
hilltops  as  normally  observed.  The  most  significant  error 
was  found  at  the  base  of  the  terrain  slope  where  winaspeeds 
were  underestimated  as  the  hydros*  :  wave  trough  was  moved 
away  from  the  mountains.  Reconunei  .  cec  improvements  include 
the  use  of  a  hydraulic  jump-like  f.  field  in  place  of  tne 
hydrostatic  wave  perturbation  in  the  model  flow  surfaces, 
and  a  more  realistic  treatment  of  the  neutral  boundary 
layer . 


f  'TViQ.-i r*  in*1* 

~T  -1  J  ‘  ”  ‘  0,1  • 

8v 


m 


|  D I  s  v 


ACKNOWLEDGEMENTS 


The  author  wishes  to  extend  sincere  thanks  and  best 
wishes  to  Peter  Lester  for  his  invaluable  insight  and  dedi¬ 
cation  in  obtaining  the  necessary  data  and  literature  which 
greatly  facilitated  this  project.  Special  thanks  also  to 
Francis  Ludwig  and  Roy  Endiich  who  made  the  tasK  of  under¬ 
standing  the  inner  workings  of  COMPLEX  much  easier.  In 
addition  the  author  would  like  to  thank  Alison  Bridger  and 
Jack  Mclodancf  for  lending  their  superior  computer  skills 
and  expertise.  Development  of  the  computer  graphics  would 
have  been  impossible  without  their  assistance.  A  very 
special  thanks  also  to  Leslie  Phillips  who  made  this  re¬ 
search  possible  by  compiling  and  digitizing  all  of  the  sur¬ 
face  wind  data  for  the  Chinook  Project.  Sincere  thanks  and 
best  wishes  to  Dan  Tredo  who  first  adapted  COMPLEX  to  the 
IBM  PC-AT  and  provided  the  author  with  invaluable  assis¬ 
tance.  Finally,  the  author  wishes  to  extend  a  warm  thank 
you  to  Donna  Hurth  for  lending  her  resourceful  and  superior 
secretarial  talents. 

The  Chinook  Mountain  Wave  Project  and  this  research 
were  partially  supported  by  a  grant  from  the  Murdock  Chari¬ 
table  Trust  of  Vancouver,  Washington. 


TABLE  OF  CONTENTS 


Page 


ABSTRACT .  Hi 

ACKNOWLEDGEMENTS  .  iv 

LIST  OF  FIGURES  .  v*  i 

LIST  OF  TABLES .  xi 

Chapter 

1 .  INTRODUCTION  .  1 

2.  LITERATURE  REVIEW  .  5 

a.  Downslope  Winds .  5 

(1)  Climatology  and  Observed 

Structure  .  5 

(2)  Causes  of  Strong  Downslope 

Windstorms .  8 

(3)  Role  of  Mesoscale  Waves  in 

Alberta  Strong  Downslope 
Windstorms .  12 


b.  Mass-Consistent  Wind  Models .  15 

3.  PROCEDURES .  23 


a.  Modification  of  COMPLEX  for  Strong 

Downslope  Windstorms  .  23 

b.  Case  and  Data  Selection .  29 

c.  General  Meteorological  Conditions.  .  31 

(1)  3  November  1975  .  31 

(2)  27  March  1986  .  35 

d.  Application  of  COMPLEX  to  SDW  Cases.  38 

(1)  Computational  Grid,  Terrain, 

and  Station  Locations  .  38 

(2)  Model  Sensitivity  Tests  ....  42 


v 


TABLE  OF  CONTENTS  (continued) 


Page 


Chapter 

4.  RESULTS .  51 

a.  Sensitivity  Tests .  51 

(1)  Effect  of  Boundary  Layer 

Thickness .  51 

'2)  Effect  of  Lee  Wave  Length  ...  56 

b.  Case  Studies .  67 

(1)  3  November  1975  .  67 

(2)  27  March  1986  .  73 

(3)  Case  Comparison  and  Discussion  .  78 

5.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 

FOR  FUTURE  RESEARCH .  8  3 

REFERENCES .  87 


LIST  OF  FIGURES 


Figure  Page 

1.  A  'volume  element*  in  the  grid  mesh 

(after  Endlich  et  al .  ,  1982) .  20 


2.  (a)  Potential  temperature  analysis  from  Boulder 

windstorm  of  11  January  1972  (after  Lilly  and 
Zipser,  1972).  Potential  temperature  isopleths 
are  drawn  at  2.5  K  intervals.  The  analysis 
extends  from  the  Continental  Divide  (CD)  through 
Boulder  (Bou) ,  Jefferson  County  Airport  (JCA) , 
and  Denver  (Den) .  (b)  Schematic  cross  section 

of  model  flow  surfaces  illustrating  terrain 
height  and  layer  thickness  parameters.  The  up¬ 
stream  tilt  of  the  hydrostatic  wave  through 

axis  is  represented  by  a  thin  dashed  line.  .  .  26 

3.  Topographical  chart  of  southwestern  Alberta 

(after  Lester,  1978a).  The  study  area  is 
indicated  by  heavy  dark  lines  while  the 
transect  line  is  indicated  by  a  dashed  line.  .  30 

4.  Synoptic  analyses  (a)  surface  analysis  at 
1400  MST  and  (b)  700  mb  analysis  at  1700  MST 
for  3  November  1975.  The  province  of  Alberta 
and  the  study  area  are  indicated  by  heavy  dark 
lines.  Both  (a)  and  (b)  incorporate  standard 
plotting  models.  In  (a;,  solid  lines  are 
isobars  at  4  mb  intervals.  The  solid  line  with 
triangular  pips  is  a  cold  front.  In  (b) ,  solid 
lines  are  contours  of  geopotential  height  at 

30  m  intervals  while  dashed  lines  are  isotherms 
at  5  C  intervals . 32 

5.  Space-time  section  of  mean  hourly  wind  speed 
(m/s)  for  3  November  1975  along  transect  line 
between  Marmot  Creek  basin  (MCB) ,  Kananaskis 
Research  Centre  (ESC) ,  Jumping  Pound  (JPd) , 
and  Calgary  International  Airport  (YYC) . 

Isotachs  are  drawn  at  5  m/s  intervals . 34 

6.  Same  as  Fig.  4,  except  for  27  March  1986..  .  .  36 


1 


Figure 

7. 

8  . 

9. 

10  . 

11 . 

12. 

13. 

14  . 


LIST  OF  FIGURES  (continued) 


Space-time  section  of  mean  hourly  wind  speed 
(m/s)  for  27  March  1986  along  transect  line 
between  Fortress  Mountain  (FtM) ,  Kananaskis 
Research  Centre  (ESC) ,  Jumping  Pound  (JPd) , 
and  Calgary  International  Airport  ( YYC ) . 

Isotachs  are  drawn  at  5  m/s  intervals . 

Three-dimensional  perspective  of  model  ter¬ 
rain  in  Calgary  area.  The  view  is  toward  the 
southwest  . 

Model  grid  domain  and  surface  stations  for 

(a)  3  November  1975  case  and  (b)  27  March  1986 

case.  Dotted  lines  are  contours  of  terrain 
height  ASL  drawn  at  100  m  intervals . 

Observed  surface  winds  for  3  November  1975 
at  1300  MST.  Stations  without  wind  vectors 
indicate  data  were  not  available  for  that 
particular  hour . 

Same  as  Fig.  10,  except  for  (a)  1400  MST  and 

(b)  1500  MST . 

Observed  surface  winds  for  27  March  1986 
at  1300  MST.  Stations  without  wind  vectors 
indicate  data  was  not  available  for  that 
particular  hour . 

Same  as  Fig.  12,  except  for  (a)  1400  MST  and 
(b)  1600  MST . 

East-west  vertical  cross  section  through 
Jumping  Pound  (JPd)  and  Calgary  International 
Airport  (YYC)  of  model  flow  surfaces  for 
(a)  Sensitivity  test  #1  and  (b)  Sensitivity 
test  #2.  Lines  from  bottom  to  top  are  the 
terrain  surface  and  the  five  upper  model  flow 
surfaces.  Note:  Sensitivity  test  #3  produced 
optimum  results.  The  corresponding  cross  sec¬ 
tion  is  presented  just  prior  to  the  case  studies 
(Fig.  22) . 


v  i  i  i 


Page 

37 

39 

40 

43 

44 

45 

46 


52 


LIST  OF  FIGURES  (continued) 


Figure  Page 

15.  Model  anemometer  height  (a)  wind  vectors  and 
fb)  isotach  analysis  for  Sensitivity  test  #1. 

In  (b) ,  isotachs  are  drawn  at  5  m/s  intervals 
and  winds  exceeding  15  m/s  are  indicated  by 
cross  hatching.  The  wave  trough  axis  is  in¬ 


dicated  by  a  dashed  line .  53 

16.  Same  as  Fig.  15,  except  for  Sensitivity 

test  #2 .  54 


East-west  variation  of  model  anemometer 
height  wind  speed  through  Jumping  Pound  (JPd) , 
Springbank  Airport  (Spr),  and  Calgary 
International  Airport  (YYC)  for  Sensitivity 
tests  #1,  #2,  and  #3.  Observed  wind  speeds 
are  denoted  by  an  asterisk  while  dotted,  solid, 
and  dashed,  lines  represent  tests  #1,  #2,  and 


#3,  respectively .  57 

18.  Same  as  Fig.  14,  except  for  (a)  Sensitivity 

cest  #4  and  (b)  Sensitivity  test  #5 .  58 


19.  Same  as  Fig.  15,  except  for  Sensitivity 

test  #4 .  59 

20.  Same  as  Fig.  15,  except  for  Sensitivity 

test  #5 .  60 

21.  Same  as  Fig.  17,  except  for  Sensitivity 
tests  #3,  #4,  and  #0.  Observed  wind  speeds 
are  denoted  by  an  asterisk  while  dashed, 
solid,  and  dotted  lines  represent  tests  if  3, 

#4,  and  #5,  respectively .  63 

22.  East-west  vertical  cross  section  through 
Jumping  Pound  (JPd)  and  Calgary  International 
Airport  (YYC)  of  model  flow  surfaces  for 
Case  Run  #1  (Sensitivity  test  #3) .  Lines  from 
bottom  to  top  are  the  terrain  surface  and  the 
five  upper  model  flow  surfaces 


65 


Figure 

23  . 

24  . 

25. 

26  . 

4.  '  . 

28  . 

29  . 

30  . 


LIST  OF  FIGURES  (continued) 


Page 


Same  as  Fig.  22,  except  for  (a)  Case  Run  #2 

and  (b)  Case  Run  #3 .  66 

Model  anemometer  height  (a)  wind  vectors  and 
(b)  isotach  analysis  for  3  November  1975  at 
1300  MST  (Case  Run  #1).  In  (b) ,  isotachs  are 


drawn  at  5  m/s  intervals .  68 

Same  as  Fig.  24,  except  for  1400  MST  (Case 

Run  #2) .  69 

Same  as  F  g.  24,  except  for  1500  MST  (Case 

Run  #3) .  70 


Model  anemometer  height  (a)  wind  vectors  and 
(b)  isotach  analysis  for  27  March  1986  at 
1300  MST  (Case  Run  #1).  In  (b) ,  isotachs  are 


drawn  at  5  m/s  intervals .  74 

Same  as  Fig.  27,  except  for  1400  MST  (Case 

Run  #2) .  75 

Same  as  Fig.  27,  except  for  1600  MST  (Case 

Run  #3) .  76 


East-west  variation  of  model  anemometer 
height  wind  speed  through  Jumping  Pound  (JPd) , 
Springbank  Airport  (Spr) ,  and  Calgary  Inter¬ 
national  Airport  (YYC)  for  modified  (dashed 
line)  and  unmodified  (solid  line)  flow  surface 
configurations  upstream  of  wave  troughs  located 
at  (a)  Btw  and  (b)  YYC.  Observed  wind  speeds 
are  denoted  by  an  asterisk  (see  text  for 
discussion) .  81 


x 


LIST  OF  TABLES 


Table  Page 

1.  Surface  wind  stations  for  (a)  3  November  1975 

and  (b)  27  March  1986  cases . 41 

2.  Observed  surface  and  700  mb  geostrophic  wind 

speeds  and  directions  for  (a)  3  November  1975 
and  (b)  27  March  1986  cases . 47 


3.  Mountain  wave  flow  parameters  assumed  for  all 
model  runs.  Values  were  derived  from  the 
potential  temperature  analysis  of  11  January 
1972  Boulder  windstorm  (Lilly  and  Zipser, 

1972) .  50 

4.  Mountain  wave  flow  parameters  specified  in 

model  sensitivity  tests . 50 

5.  Mountain  wave  parameters  specified  for 

3  November  1975  and  27  March  1986  case  studies. 64 


XT 


Chapter  1 


INTRODUCTION 


a .  Background 

It  is  well  known  that  the  flow  of  stable  air  over  moun¬ 
tain  barriers  produces  significant  disturbances  through  deep 
layers  downwind  of  the  mountains.  On  the  mesoscale,  moun¬ 
tain-induced  disturbances  include  lee  waves,  rotor  circula¬ 
tions  ,  and  dowr.slope  windstorms.  Over  the  past  5C  years, 
many  field  studies,  theoretical  investigations,  and  numeri¬ 
cal  modeling  studies  have  contributed  significantly  to  the 
understanding  of  these  pnencmena  (e.g.,  Lyra,  1943;  Queney, 

1 9 4 S ;  Scorer,  1949;  Long,  1954;  Alaka,  1960  ;  Nichclls,  1973  ; 
Kuettner  and  O'Neill,  1981;  Lilly  et  al . ,  19S2;  Smith, 

19  3  5) 

Dowr.slope  windstorms,  known  in  many  areas  of  the  world 
by  local  names,  are  either  warm  (e.g.,  Alpine  Foehn,  Rocky 
Mountain  Chinook,  Argentine  Ccnda)  or  cold  (e.g.,  French 
Mistral,  Yugoslavian  Bora)  dry  and  gusty  winds  from  the 
mountains  (Huschke,  1959).  The  importance  of  downslope 
windstorms  lies  in  their  role  in  the  global  angular  momentum 
balance  (Bretherton,  1969;  Lilly,  1972)  and  in  potentially 
serious  lr  effects.  For  example,  downslope  windstorms 
occasion.  •  develop  to  hurricane-strength  (_>  35  m  s  ±)  . 
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These  so  called  strong  dcwnslope  windstorms  (SDW)  are 
capable  of  producing  widespread  damage  and  severe  aircraft 
turbulence  (Lilly  and  Zipser,  1972;  Fingerhut  and  Lester, 
1973;  Lilly,  1978;  Lester,  1978a;  Hoinka,  1985;  Lester  and 
Bach ,  1986 ) . 

In  southern  Alberta,  Chinooks  are  a  common  occurrence 
along  the  eastern  slopes  of  the  Canadian  Rocky  Mountains 
(Longley,  1967;  Lester,  1975,  1976a,  1976b;  Giusti,  1987). 
The  warm  Chinook  winds  act  to  ameliorate  an  otherwise  severe 
winter  climate.  Population  and  agriculture  distributions 
clearly  reflect  the  presence  of  the  Chinook  belt  and  it 
follows  that  accurate  prediction  of  chinook  characteristics 
is  an  important  task  for  forecasters  in  that  area. 

The  state  of  the  art  of  chinook  forecasting  in  Alberta 
is  similar  to  the  status  of  prediction  capabilities  else¬ 
where  iScheet:  et  al .  ,  1976;  Sangster,  1977;  Brown,  1986). 
Synoptic  scale  flow  patterns  related  to  chinook  occurrence 
are  well  documented  and  can  be  predicted  with  a  high  degree 
of  accuracy  (Lester,  1976a;  Drews,  1982),  but  mesoscale 
problems  still  exist.  For  example,  the  observational  net¬ 
work  in  southern  Alberta  cannot  resoive  mesoscale  character¬ 
istics  of  chinook  airflow  (Lester,  1976a).  This  is  a 
serious  problem  since  the  worst  chinook  windstorms  are 
driven  by  mesoscale  processes  (Lilly  and  Zipser,  1972; 

Lilly,  1978;  Lester,  1978a).  Additionally,  the  theory  of 
downslope  windstorms  is  incomplete  (Durran,  1986a)  and 


current  numerical  models  are  either  oversimplified  or 
require  computer  resources  which  regional  forecast  offices 
do  not  possess  (Peltier  and  Clark,  1979;  Durran,  1986a; 
Lester,  1987).  As  a  result,  precise  chinook  windstorm 
details,  such  as  the  timing  of  the  onset  and  the  loca¬ 
tion/magnitude  of  the  strongest  winds,  cannot  be  accurately 
predicted.  This  problem  is  magnified  in  the  case  of  SDWs. 

A  partial  solution  to  the  SDW  forecast  problem  has  been 
found  in  some  areas  through  the  use  of  special 
’ mesonetworks '  of  observation  stations.  Subsequent  analyses 
of  data  from  such  networks  (e.g.,  Lilly  and  Zipser,  1972; 
Brinkmann,  1974;  Lilly,  1978;  Lester,  1978a;  Bedard  et  al . , 
1981;  Mathews  et  al . ,  1984)  have  suggested  that  SDWs  often 
behave  coherently  on  the  mesoscale,  i.e.,  onsets  and  move¬ 
ments  cf  the  windstorms  may  be  extrapolated  over  time  and 
space  scales  of  a  few  hours  and  several  tens  of  kilometers. 
These  results  indicate  that  regular  mesoscale  'nowcasts'  and 
short  range  forecasts  are  possible  given  detailed  observa¬ 
tions  of  mesoscale  flow. 

One  difficulty  with  this  solution  in  Alberta  is  that 
network  costs  would  be  prohibitively  high  because  of  the 
large  area  of  the  chinook  belt  and  the  complexity  of  the 
terrain.  It  is  proposed  here  that  this  difficulty  can  be 
overcome  by  combining  a  moderately  dense,  but  affordable 
station  network  with  an  accurate,  fast,  and  economical 
interpolation  scheme  which  takes  into  account  the  complexity 
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of  the  terrain.  Recent  advances  in  objective  wind  analysis 
techniques  for  complex  terrain  indicate  that  this  proposal 
is  feasible  (e.g.,  Enalich,  1967;  Sherman,  1978;  Bhumralkar 
et  al.,  1978;  Endlich  and  Lee,  1983;  Ludwig  et  al.,  1985; 
Tredo,  1987).  In  addition  to  these  developments,  the 
Alberta  Chinook  Project;  (Lester,  1985  ;  Lester  and  Phillips, 
1987)  has  produced  a  series  of  windstorm  data  sets  which  are 
suitable  for  mesoscale  analysis.  Thus,  the  objective  of 
this  study  is  to  investigate  the  feasibility  of  analyzing 
severe  chinook  windstorms  in  Alberta  with  a  diagnostic, 
boundary  layer  wind  model 

The  investigation  is  carried  out  by  (i)  modifying  a 
state-of-the-art,  diagnostic,  mass-consistent  wind  model 
(COMPLEX)  to  take  into  account  the  influence  of  mountain  lee 
waves,  (ii)  testing  the  sensitivity  of  the  model  to  those 
modifications,  and  (iii)  applying  the  model  to  two  SOW 
cases.  In  Chapter  2,  observations  and  theory  of  downslope 
winds  and  mass  consistent  wind  models  for  complex  terrain 
are  reviewed  to  provide  the  background  for  the  modification 
of  COMPLEX  for  mountain  wave  effects.  A  description  of  the 
procedures  used  to  modify  COMPLEX,  the  case  data,  sensi¬ 
tivity  tests,  and  experimental  procedures  are  presented  in 
Chapter  3.  Results  are  discussed  in  Chapter  4,  and  a  brief 
summary  of  results,  conclusions,  and  recommendations  for 
future  research  are  presented  in  Chapter  5. 
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Chapter  2 


LITERATURE  REVIEW 

a .  Downslope  Winds 

( l )  Climatology  and  Observed  Structure 
Strong  downslope  windstorms  (SOW)  are  a  subset  of  oro¬ 
graphic  flows  in  which  hurricane  force  winds  are  produced 
along  the  lee  slopes  of  mountains  (Brinkmann,  1974;  Smith, 
1^85; .  Foehn  and  Bora  represent  the  two  general  types  of 
downslope  winds,  which  are  traditionally  classified  accord¬ 
ing  to  temperature  effects  [Atkinson,  1981;  Durran,  1986a). 
While  SDW  may  occur  with  both  Bora  and  Foehn,  the  former 
cause  temperatures  to  decrease  while  the  latter  produce  tem¬ 
perature  increases. 

Chinook  winds  are  of  Foehn  type  (Brinkmann,  1973). 

They  occur  along  the  entire  length  of  the  Rocky  Mountain 
chain  and  have  been  studied  in  detail  in  southern  Alberta 
(Longley,  1967;  Brinkmann  and  Ashwell,  1968;  Holmes  and 
Hage,  1971;  Kellie,  1972;  Lester,  1975,  1976a;  Giusti ,  1987) 
and  in  Colorado  (Ives,  1950;  Julian  and  Julian,  1969;  Miller 
et  al . ,  1974;  Brinkmann,  1974;  Scheetz  et  al . ,  1976). 

The  Chinook  of  southern  Alberta  is  well-defined  because 
the  Arctic  front  with  its  extreme  temperature  contrasts 
often  lies  in  that  area  (Mathews  et  al.,  1986).  As  storms 
move  across  the  area,  the  alternation  between  cold  arctic 
air  and  warm  chinook  conditions  is  dramatic  and  frequent. 
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Lester  (1975,  1976b)  has  quantified  the  traditional  descrip¬ 
tion  of  the  Chinook  (i.e.,  a  warm,  dry,  gusty  wind  from  the 
mountains  associated  with  rapid  temperature  changes)  to 
derive  a  Chinook  climatology  for  Calgary,  located  in  south¬ 
western  Alberta.  From  five  years  of  data  (1971-1976),  he 
found  that  Chinooks  occurred  on  an  average  of  13.5  days  each 
month  during  October-Mar ch .  A  'typical'  (median)  Alberta 
Chinook  had  a  duration  of  8  h,  minimum  relative  humidity 
near  40%,  wind  direction  WSW-WNW,  and  hourly  averaged  wind 
speed  of  8  m  s-1  with  maximum  gusts  to  17.5  m  s~^. 

Lester  (1976b)  showed  that,  for  Calgary,  the  frequency 
of  'extreme'  Chinook  windstorms  (i.e.,  those  with  gusts  >_  25 
m  s~~)  was  8%  between  October  and  March.  For  the  same 
months,  but  for  1978-1984  at  a  location  in  the  foothills 
west  of  Calgary,  Giusti  (1987)  determined  that  13%  of 
Chinook  windstorms  could  be  classified  as  'extreme1.  He 
found  that  windstorms  occurred  more  frequently  during  the 
afternoon,  although  this  tendency  was  less  distinct  for  the 
strongest  windstorms.  Also,  windstorm  durations  were  less 
than  four  hours  for  more  than  half  the  cases. 

Individual  Chinook  windstorms  are  characterized  by 
their  intermittency  in  space  and  time,  i.e.,  variations  in 
movement  and  gradients  in  intensity.  Records  of  wind  speed 
from  severe  Chinook  windstorms  (Lilly  and  Zipser,  1972; 
Brinkmann,  1974;  Lester,  1978a;  Mathews  et  al.,1986)  show 
that  space  and  time  scales  of  the  chinook  winds,  SDW,  and 
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wind  gusts  span  a  wide  range.  On  the  largest  scales  are  the 
Chinook  winds  (i.e.,  warm,  dry,  winds  from  the  mountains 
with  speeds  of  several  m  s-^) .  The  area  of  influence  can 
extend  as  far  as  350  km  away  from  the  Continental  Divide 
(Lor.gley,  1967;  Mathews  et  al.,  1984)  and  over  that  distance 
along  the  divide.  These  general  Chinook  conditions  may  last 
a  few  days  (Lester,  1976b).  In  contrast,  severe  cninook 
windstorms  which  include  SDW  have  a  reduced  area  of  influ¬ 
ence  and  a  shorter  duration.  In  southern  Alberta,  hourly 
averaged  wind  speeds  exceeding  20  m  s-^  typically  occur  in  a 
rone  40-50  km  wide  near  the  base  of  the  terrain  slope  and 
last  3-4  hours  (e.g.,  Lester,  1978a).  On  the  smallest 
scales  are  windy  periods  of  tens  of  minutes  and  individual 
gusts  which  may  last  on  the  order  of  seconds. 

The  observed  movements  of  SDWs  are  quite  variable 
(Lester  and  Phillips,  1987).  In  southern  Alberta  and  in 
Boulder,  windstorm  locations  have  been  observed  to  fluctuate 
up  and  down  the  slopes  and  to  move  steadily  out  into  the 
plains  toward  the  east  (Brinkmann.  1974;  Lester  1978a).  In 
both  regions,  gradients  in  windstorm  intensity  are  highest 
at  the  base  of  the  terrain  slope.  Lester  (1976a)  reported 
that  the  onset  of  the  chinook  in  southern  Alberta  is  often 
marked  by  wind  speed  maxima  immediately  to  the  lee  of  the 
mountains.  A  similar  observation  was  reported  for  Boulder 
by  Brinkmann  (1974)  from  her  analysis  of  twenty  cases  in 
which  gusts  exceeded  33  m  s-^.  A  dominant  feature  of  the 
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Boulder  surface  winds  was  a  large  peak  in  wind  speed  at  the 
foot  of  the  mountains  with  rapid  decreases  (50%  over  a  dis¬ 
tance  of  10  km)  to  the  east  and  west. 

( 2 )  Causes  of  Strong  Downslope  Windstorms 

Severe  Chinook  windstorms  are  the  result  of  the  inter¬ 
action  of  circulations  with  a  wide  range  of  horizontal 
scales  (Lester,  1975,  1978a;  Lilly,  1978).  For  example,  in 
Aioerta,  chinock  windstorms  occur  when  the  planetary  wave 
pattern  produces  strong,  nearly  zonal  tropospheric  flow  of 
anomalously  warm  air  across  the  Canadian  Rocky  Mountains. 
This  cross-mountain  flow  is  intensified  by  embedded,  rapidly 
moving,  barociinic  wave  disturbances.  These  conditions  are 
manifested  by  an  intense  surface  cyclone  in  the  Gulf  of 
Alaska  and  a  ridge  of  high  pressure  building  in  the  state  of 
Washington.  Low  pressures  dominate  in  northern  British 
Columbia  (B-C.)  and  along  the  lee  of  the  Rockies  in  Alberta 
(Lester,  1976b).  As  a  frontal  system  crosses  B.C.  and 
approaches  the  Continental  Divide,  pressure  falls  along  the 
lee  slopes  and  rises  in  southeastern  B.C.  These  events  lead 
to  the  formation  of  a  ridge  of  sea  level  pressure  over 
southern  B.C.  (Brinkmann,  1969),  a  common  synoptic  feature 
of  the  Chinook  known  locally  as  the  "Foehn  nose". 

On  the  mesoscale,  the  causes  of  downslope  windstorms 
are  tightly  linked  to  large-amplitude  mountain  lee  waves 
forced  by  the  previously  discussed  synoptic  flow  conditions 
fe.g.,  Long,  1954;  Houghton  and  Kasahara,  1968;  Klemp  and 
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Lilly,  1975,  1978;  Clark  and  Peltier,  1977;  Peltier  and 
Clark,  1979;  Blumen  and  Hartsough,  1985;  Durran,  1986a, 
1986b) .  Mountain  lee  waves  are  stationary  gravity  waves 
which  form  downwind  when  the  stable  airflow  is  displaced 
vertically  by  an  orographic  barrier  (Durran,  1986a). 

Over  the  past  30  years,  three  theories  have  been 
advanced  to  explain  the  role  of  lee  wave  phenomena  in  SDW 
development.  These  are  (i)  hydraulic  jumps,  (ii)  partial 
reflection  of  upward  propagating  waves,  and  (iii)  trapping 
of  energy  by  self-induced  critical  layers. 

According  to  the  hydraulic  jump  theory  (Long,  1954; 
Houghton  and  Kasahara,  1968),  strong  winds  develop  along  the 
lee  slope  as  the  fluid  undergoes  a  transition  from  upstream 
subcritical  flow  to  supercritical  flow  over  the  mountain 
(Durran,  1986a).  The  strength  of  this  theory  is  that  the 
nonlinearity  introduced  by  finite  terrain  barriers  is  ex¬ 
plicitly  accounted  for.  However,  hydraulic  jump  theory 
suffers  from  the  unrealistic  requirement  that  the  fluid  top 
must  be  bounded  by  a  rigid  lid  or  free  surface.  As  a  re¬ 
sult,  the  vertical  propagation  of  energy  which  is  observed 
to  occur  in  large  amplitude  mountain  waves  cannot  be  in¬ 
cluded  in  this  theory. 

The  partial  reflection  of  vertically  propagating  waves 
as  a  factor  in  SDW  development  was  investigated  by  Klemp  and 
Lilly  (1975,  1978)  using  a  three-layer,  linear,  two-dimen¬ 
sional,  hydrostatic  model  with  bell-shaped  terrain  (Durran, 
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1986b) .  Upward  propagating  waves  are  influenced  by  the  wind 
and  stability  profiles  as  defined  by  the  Scorer  parameter 
( $2  )  r  where 

J?2  =  (N/U  )2  -  1/U  (d2U/'dz2 )  ,  (1) 

N2  =  (g/60) (de/dz)  ,  (2) 

N  is  Brunt-Vaisala  frequency,  U  is  mean  horizontal  wind 

speed,  g  is  the  gravitational  constant,  ©0  is  a  reference 

potential  temperature,  and  (de/'dz)  is  the  mean  vertical 

potential  temperature  gradient.  If  the  wave  encounters  a 
I  '■* 

region  m  which  a*  changes  rapidly  with  height,  part  of  the 
wave  energy  may  be  reflected  downward.  As  a  result,  the  wave 
amplitude  below  the  reflecting  layer  is  determined  by  the 
superposition  cf  the  upward  and  downward  propagating  waves. 
For  the  special  case  of  a  three-layer  atmosphere  flowing 
over  sinusoidal  topography,  Kiemp  and  Lilly  found  that  the 
optimal  superposition  occurs  when  the  vertical  dimension  of 
each  lower  layer  is  one-fourth  of  the  vertical  wavelength. 
Applying  a  multi-layered  numerical  model  to  the  Boulder 
region,  they  achieved  reasonable  success  in  predicting 
downslope  winds  from  upstream  soundings.  Their  multiple 
layer  approach  has  the  advantage  that  with  a  sufficient 
number  of  layers,  the  mean  state  can  be  configured  for  a 
close  match  with  an  actual  sounding.  Blumen  and  Hartsough 
(1985)  improved  the  Kiemp  and  Lilly  approach  with  a  continu¬ 
ous  model  for  partial  reflection.  However,  it  must  be  kept 
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in  mind  that  because  the  basic  Klemp-Lilly  model  is  linear, 
caution  must  be  exercised  in  its  application  to  large  ampli¬ 
tude  waves  where  the  perturbation  fields  are  comparable  in 
magnitude  to  the  mean  flow. 

Nonlinear  and  non-hvdrostatic  effects  from  large-ampli¬ 
tude  mountain  waves  on  SDW  development  were  investigated  by 
Clark  and  Peltier  (1977)  and  Peltier  and  Clark  (1979)  using 
a  two-dimensional  model  with  bell-shaped  terrain  (Durran, 
1986b).  They  concluded  that  if  the  mountain  wave  amplifies 
and  overturns  (breaks)  just  above  the  tropopause,  the  wave 
energy  becomes  trapped  between  the  ground  and  the  level  of 
breaking.  The  overturning  region  is  interpreted  as  a  self- 
induced  critical  layer,  and  is  characterized  by  the  produc¬ 
tion  of  locally  superadiabatic  regions,  strong  mixing,  and  a 
local  wind  reversal.  The  trapping  of  the  wave  energy  pro¬ 
duces  a  large  amplitude  wave  which  efficiently  transfers 
large  amounts  of  momentum  to  the  surface,  producing  a  SDW. 
Their  approach  explicitly  accounts  for  nonlinear  and  non¬ 
hydrostatic  effects,  but  it  has  the  disadvantage  that  the 
interpretation  of  the  complex  interaction  of  several  pro¬ 
cesses  is  sometimes  difficult.  Although  controversy  still 
exists  about  the  role  of  self-induced  critical  layers  on 
wave  amplification  (Lilly  and  Klemp,  1980;  Peltier  and 
Clark,  1980),  it  is  generally  accepted  that  non-linear 
effects  are  important  (Durran,  1986a) . 
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(3)  Role  of  Mesoscale  Waves  in  Alberta  Strong 
Pownslope  Windstorms 

Although  no  comprehensive  measurements  have  been  made 
of  the  three-dimensional  mesoscale  structure  of  SDWs  in 
Alberta,  much  indirect  evidence  has  been  gathered  in  support 
of  a  significant  role  for  mountain  lee  waves  in  SDW  produc¬ 
tion  in  that  area.  From  a  12  year  record  of  lenticular  and 
rotor  cloud  ooservations  at  Pincher  Creek  in  southwestern 
Alberta,  Lester  <19C8b)  deduced  a  very  high  frequency  of  lee 
waves.  He  found  that  lee  wave  clouds  occur  on  an  average  of 
141  days  annually  and  on  50%  or  more  of  the  days  during  the 
most,  active  months  (October  and  February). 

The  presence  of  50-100  km  lee  waves  (i.e.,  those  which 
are  likely  to  be  hydrostatic)  during  Chinook  situations  was 
first  documented  by  Holmes  and  Hage  (1971).  Using  instru¬ 
mented  aircraft,  they  found  that  large  snowmelt  areas  were 
produced  in  the  troughs  of  waves  that  had  carried  warm  air 
downward,  displacing  very  cold  air  at  the  surface.  In  a 
particular  case  (13  January  1963),  the  wavelengths  ranged 
from  64-67  km  and  the  wave  amplitude  was  estimated  at  500- 
7  00  m . 

Other  evidence  of  long  lee  waves  in  the  area  was  given 
by  Kellie  (1972)  who  applied  spectrum  analysis  to  data  from 
two  instrumented  towers  in  Calgary.  He  found  eastward  pro¬ 
gressing  waves  with  lengths  of  29  to  118  km. 
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Lester  (1S75,  1976a,  1976c)  and  Lester  and  McPherson 
(1977)  have  used  aircraft,  satellite,  and  surface  observa¬ 
tions  used  to  show  that  the  "chinook  arch,"  a  distinctive 
lenticular  cloud  that  often  extends  hundreds  of  kilometers 
along  the  lee  of  the  mountains,  lies  in  the  crest  of  a 
gravity  wave  with  a  length  of  5G-100  km.  Furthermore,  they 
found  that  the  waves  often  propogated  eastward  from  the 
mountains  with  a  phase  speed  less  than  the  windspeed. 
Occasionally ,  observations  of  multiple  chinook  arches  sug¬ 
gested  the  presence  cf  at  least  two  long  waves  downwind  of 


the  mountains. 

Lester  (197 

6d)  documented 

the  high 

frequency  of  occur 

r en c e  of  chinook 

arches  land. 

therefore , 

very  long  lee 

waves)  observed  from  Calgary  during  16  months  in  1974-1976. 
Arches  were  observed  on  132  days  with  frequencies  as  high  as 
15  days  per  month  in  December  and  January.  He  also  found 
that  the  relationship  between  arch  occurrence  and  gusty 
chinook  winds  was  very  strong  m  the  mountains,  but 
decreased  significantly  to  the  east.  For  example,  at 
F.ananaskis  Forestry  Experiment  Station,  about  10  km  east  of 
the  Continental  Divide,  92%  of  the  arch  clouds  (25  cases) 
occurred  with  gusty  chinook  winds,  while  at  Calgary,  the 
percentage  decreased  to  48%.  At  this  time,  no  clear  rela¬ 
tionship  has  been  established  between  the  appearance  of  the 
chinook  arch  and  the  occurrence  of  SDW. 
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In  a  recent  climatological  study,  Giusci  (1987)  exam¬ 
ined  soundings  taken  immediately  upstream  of  the  Continental 
Divide  just  prior  to  Chinook  windstorms.  He  stratified 
static  stability  profiles  by  observed  Chinook  wind  speed 
categories  for  several  years  of  data  and  found  that  the 
strength  of  the  wind  storm,  increased  with  an  increase  of 
static  stability  at  mountain  top.  A  similar  structure  has 
been  associated  with  Colorado  windstorms  ( Br inkm.ann  ,  IS-  3 
and  agrees  with  theoretical  stability  structures  that  Klemp 
and  hilly  (1975)  have  found  to  be  conducive  to  SDW. 

The  evidence  for  the  importance  of  nor.-hydrostatic 
waves  and  non-linear  effects  in  the  production  of  chinook 
windstorms  ir.  Alberta  is  not  as  extensive  as  that  for  the 
hydrostatic  lee  wave.  A  recent  paper  by  Lester  and  Bach 
■1 986)  examined  an  extreme  clear  air  turbulence  incident 
wr.ich  occurred  near  the  tropcpause  during  an  SDW  event 
approximately  50  km  west  of  Calgary.  Satellite  data  and 
derived  ait  craft  records  of  vertical  velocity  and  potential 
temperature  showed  that  both  long  (50-100  km)  hydrostatic 
lee  wave,  and  shorter  (10-20  km)  non-hydrostatic  mountain 
lee  waves  we’-e  present.  Furthermore,  an  extensive  turbulent 
layer  which  existed  near  the  tropopause  was  suggestive  of 
the  critical  layer  necessary  for  wave  energy  trapping,  as 
hypothesized  by  Peltier  and  Clark  (1979). 
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b. 


Mass-Consistent  Wind  Models 


Even  with  recent  advances  in  downslope  wind  theory, 
analysis  of  surface  wind  fields  associated  with  SDWs  remains 
a  difficult  problem  (e.g.,  Blumen  and  Hartsough,  1985; 

Smith,  1985;  Durran,  1986a).  The  numerical  models  developed 
by  Klemp  and  Lilly  (1975),  Peltier  and  Clark  (1979),  and 
Durran  (19S6a)  simulate  the  dynamics  of  mountain  airflow, 
and  can  produce  mountain  wave  patterns  consistent  with 
observations.  However,  their  use  of  smoothed  topography 
limits  direct  application  of  those  models  to  boundary  layer 
wind  flow  because  smaller  scale  flow  perturbations  due  to 
complex  terrain  are  not  included.  As  noted  in  Chapter  1, 
the  current  investigation  will  attempt  to  overcome  the 
latter  problem  by  the  adaptation  of  a  diagnostic  wind  model 
to  yield  physically  reasonable  hourly  objective  surface  wind 
analyses  during  SDW  events. 

Mass-consistent  wind  models  (e.g.,  Sherman,  1978; 
Bhumralkar  et  al . ,  1978)  are  objective  boundary  layer  wind 
analysis  schemes  which  are  applicable  to  regions  of  complex 
terrain  and  require  only  modest,  readily  available  computer 
resources  (Tredo,  1987).  In  the  most  sophisticated  of  these 
models,  wind  observations  are  interpolated  to  a  three-dimen¬ 
sional  grid  and  then  adjusted  so  that  conservation  of  mass 
is  satisfied.  The  accuracy  of  the  final  wind  field  is 
dependent  upon  the  number  and  quality  of  wind  observations, 
and  on  the  representation  of  physical  processes  by  the  model 
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inputs  (Mass  and  Dempsey,  1985;  Tredo,  1987). 

One  of  the  earliest  diagnostic  wind  models  was  the 
mass-adjusted,  three-dimensional  wind  flow  model  for  complex 
terrain  (MATHEW)  developed  by  Sherman  (1978)  from  the  conti¬ 
nuity  equation  in  rectangular  coordinates  (Ludwig,  1987). 

In  MATHEW,  a  variational  method  by  Sasaki  (1958,  1970)  was 
used  to  minimize  differences  between  the  initial  wind  analy¬ 
sis  derived  from  observed  values  and  the  final  mass-consis¬ 
tent  flow.  A  power  law  was  used  for  vertical  profile  esti¬ 
mates.  Practical  use  of  MATHEW  is  difficult  because  the 
rectangular  coordinates  complicate  the  problem  of  specifying 
meteorologically  satisfactory  boundary  conditions,  espe¬ 
cially  at  the  lower  boundary  (Endlich,  1967;  Ludwig,  1987). 

In  the  COMPLEX  mass-consistent  wind  model  developed  by 
Bhumralkar  et  al.  (1978),  the  sigma  (c)  system  is  incorpo¬ 
rated  to  produce  an  initially  terrain-following  coordinate 
system . 


c  =  [  z  -  h(x,y!J/[H(x,y)  -  h ( x , y ) ]  ,  (3) 

where  z  is  height  above  sea  level,  h  is  height  of  terrain 
above  sea  level,  and  H  is  the  height  of  the  boundary  layer 
top.  The  lower  boundary  condition  is  specified  by  the 
terrain  height  while  the  upper  boundary  (boundary  layer  top) 
is  given  by  the  equation; 

H(x,y)  =  Z  +  sh(x,y)  +  (1  -  s)hQ  ,  (4) 
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where  Z  is  the  average  boundary  layer  tnickness,  s  is  a 
"slope"  factor,  which  defines  the  degree  to  which  the  bound¬ 
ary  parallels  the  underlying  terrain,  and  hQ  is  tne  terrain 
height  at  a  specified  reference  point.  The  use  of  sigma  co¬ 
ordinates  allows  for  simple  formulation  of  the  upper  and 
lower  boundary  conditions  as  the  terrain  and  boundary  layer 
top  are  assumed  to  be  impermeaole  boundaries.  Observations 
cf  horizontal  wind  from  scattered  stations  are  interpolated 
to  a  horizontal  grid  and  extrapolated  vertically  to  each 
sigma  level. 

More  recently,  COMPLEX  has  been  modified  to  include  (i; 
the  use  cf  the  gecstrophic  wind  to  represent  upper  boundary 
flow  (Endlich  et  ai.,  1930,  1932,',  iii|  tne  incorporation  of 
critical  streamline  concepts  to  convert  from  a  terrain- 
following  to  a  flow-following  coordinate  system  (Endlich  and 


Let ,  . 9 o j ;  Luc wig  et 

al . ,  1935 i , 

and 

!iii)  the 

incius i on  of 

an  iterative  (rather 

than  variati 

onal 

calculus ) 

procedure  to 

eliminate  divergence 

in  the  wind 

*•  ~  o 

d  (Endlich 

,  1 9 1 "  ; 

Endlich  et  ai . ,  193C ) 

These  changes 

are  aescr 

ibed  below. 

The  gecstrophic  wind  is  applied  at  a  specified  maximum 
heignt  above  the  terrain  unless  the  boundary  layer  top  is 
below  the  mountain  tops,  so  that  the  influence  of  the  under¬ 
lying  terrain  is  limited  to  a  shallower  layer.  The  critical 
streamline  concept  (Hunt  and  Snyder,  1930)  involves  the 
balance  between  the  potential  energy  gained  from  vertical 
motions  in  a  stably  stratified  atmosphere  and  the  kinetic 
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energy  of  the  horizontal  motion  (Ludwig  et  al . ,  1985; 

Ludwig,  1987;  Tredo ,  1987).  At  the  critical  streamline 
height,  kinetic  energy  of  the  flow  equals  the  potential 
energy  gained  by  raising  air  parcels  to  the  top  of  a  terrain 
obstacle.  The  governing  relationship  for  a  stable,  adia¬ 
batic  atmosphere  in  which  vertical  displacements  are  small 
can  be  written  as 


(5; 


The  left  side  is  the  kinetic  energy  (per  unit  mass)  of  the 
flow  while  the  rignt  side  is  the  potential  energy  (per  unit 
mass;  gained  by  lifting  the  air  from  ZQ  to  the  height, 

66/5z  is  the  vertical  gradient  of  potential  temperature,  and 
T_  is  the  average  temperature  through  the  layer.  Integrat¬ 
ing  the  right  side  of  Eqn.  (5)  gives  the  highest  rise  for 
the  flow  against  the  existing  lapse  rate, 


—  t  r  r  /  _  /  rr» 


[  (S/T_)  ( 5 0 / 5 z ) n ] 


where  (  )Q  refers  to  conditions  at  the  lowest  geometric 
level.  Air  above/below  the  critical  streamline  height  pro¬ 
ceeds  over/around  the  obstacle,  provided  that  wind  increases 
with  height.  In  the  model,  the  highest  obstacles  over  which 
the  air  can  flow  are  defined.  Then  sigma  levels  are  reas¬ 
signed  from  terrain-following  to  flow-following  based  on 


18 


this  critical  streamline  concept.  Wind  speed  over  the  lower 
terrain  elevations  and  lapse  rates  at  the  various  sigma 
levels  form  the  basic  for  rh: c  definition  of  flow  surface 
shape . 

Once  the  flow-following  surfaces  are  determined,  diver¬ 
gence  is  eliminated  using  an  iterative  procedure  (Endlich, 
1967;  Endlich  et  al .  ,  1980)  in  which  horizontal  wind  compc- 


nents  on  volume 

elements 

'rig-  U 

are  sequentially 

adjusted 

to  balance  mass 

inflow  wi 

rh  outfl 

ow  (Tredo ,  1987 )  . 

n  a  c  r. 

volume  element  t 

s  bounded 

at  the 

top  and  bottom  by 

flow  sur 

faces  so  that  inflow  or  outflow  from  each  element  can  occur 
only  through  the  vertical  faces.  The  divergence  (D)  for  a 
volume  element  is  given  by 

D  =  ;5u*  ox)  +  (5v*/oy)  ,  (7) 

u  *  —  u  5  c  ,  ( S ) 

V*  =  V  X  5c  ,  (9) 

where  u  and  v  are  the  hortcontai  velocity  components  on  the 
volume  element;  6c  is  the  vertical  separation  between  two 
flow  surfaces,  and  u*  and  v*  are  the  respective  fluxes  for 
each  hortcontai  velocity  component.  For  example,  if  conver¬ 
gence  is  occurring  in  a  volume  element  due  to  greater  in¬ 
flow,  wind  components  are  adjusted  to  deer ease/ increase  wind 
component  values  on  inflow/outflow  faces  of  the  volume 
element  until  the  value  of  D  is  zero.  However,  this  adjust¬ 
ment  changes  the  divergence  value  for  adjacent  volume 
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elements,  therefore  the  procedure  is  repeated  until  the 
divergence  in  all  volume  elements  is  removed.  In  practice, 
the  iterations  are  stopped  after  a  specified  number  have 
been  completed,  resulting  in  small  divergence.  By  making 
the  wind  field  non-diver gent ,  those  boundary  conditions  tha 
are  usually  required  to  solve  Poisson  equations  for  stream 
function  and  vorticity  are  net  needed. 

The  final  step  is  to  interpolate  the  mass-adjusted  win 
components  on  flow-following  surfaces  to  anemometer  height 
(10  m)  above  the  terrain.  A  detailed  discussion  cf  the  co¬ 
ordinate  system,  initial  wind  analyses,  mass -conservation 


ad jus  trnent 

,  and 

surface 

-wind 

interpolation  procedure  is  pr e 

sen  ted  by 

T-  edc 

(193“  ;  . 

Ste 

ps  in  the  computational  prece- 

dure  (Nitz 

et  al 

.,  1935 

;  are 

as  follows: 

i;  Specify  boundary  layer  top  and  terrain-following 
(sigma,'  surfaces. 

2;  Interpolate  wind  profile  and  surface  winds  to 
sigma  surfaces. 

3/  Make  initial  (first  guess)  objective  wind  analysis 

4)  Adjust  sigma  surfaces  to  become  flow-following 
surfaces  using  the  critical  streamline  concept. 

5)  Remove  divergence  from  wind  field. 

6 1  Interpolate  mass-adjusted  wind  field  to  anemometer 
height  (10  m)  above  the  terrain  surface. 

Results  from  investigations  of  Chinook  boundary  layer 
structure  using  acoustic  remote  sensors  (SODARS)  near 
Calgary  (Mathews  and  Hicks,  1079;  Mathews  et  ai.,  1984) 
indicate  that  the  onset  and  cessation  of  the  chinook  wind¬ 
storms  is  often  strongly  influenced  by  the  interaction  of 
mountain  waves  with  the  surface  boundary  layer  (Lester, 
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1935) .  Since  COMPLEX  has  evolved  to  a  point  where  it  can 
produce  mass-consistent  wind  fields  for  varied  airflow 
structures  (e.g.,  mountain  waves),  its  application  to 
Alberta  Chinook  SDWs  is  an  appropriate  experiment  which  has 
not  been  previously  attempted. 
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Chapter  3 


PROCEDURES 

a .  Modification  of  COMPLEX  for  Strong  Downsiope  Windstorms 

Despite  the  past  successes  of  COMPLEX  in  the  analysis 
of  winds  in  mountainous  terrain,  the  initial  application  of 
that  model  to  a  Chinook  windstorm  produced  unrealistic 
surface  wind  patterns.  For  example,  the  strongest  winds 
occurred  only  over  the  mountain  peaks  with  little  enhance¬ 
ment  of  the  winds  at  the  foot  of  the  mountains  as  is  usually 
observed.  This  result  was  due  to  the  fact  that  COMPLEX  (as 
described  in  Chapter  2)  dees  not  take  into  account  the 
important  influence  of  mountain  lee  waves  on  the  boundary 
layer.  In  this  section,  the  modification  of  COMPLEX  to 
reflect  the  influence  of  lee  waves  is  described. 

As  noted  in  Chapter  2,  the  critical  streamline  concept 
(Hunt  and  Snyder,  1980)  has  been  used  in  recent  versions  of 
COMPLEX  to  make  the  transition  from  a  terrain-following  to  a 
f low- f o blowing  coordinate  system  (e.g..  Trade,  1987).  While 
this  modification  introduces  more  reasonable  physics  into 
COMPLEX,  it  does  not  incorporate  mesoscale  forcing  of  sur¬ 
face  winds  by  mountain  waves.  In  addition,  if  the  assump¬ 
tions  underlying  the  critical  streamline  concept  do  not 
hold,  then  the  use  of  flow-following  coordinates  by  COMPLEX 
requires  knowledge  of  the  three-dimensional  structure  of  the 
flew,  including  the  height  of  the  boundary  layer  top. 
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Unfortunately,  such  observations  are  not  available  for 
Chinook  SOW  cases  considered  here.  Although  the  Chinook 
Project  has  generated  several  windstorm  data  sets  over  the 
last  few  years,  these  sets  are  mainly  composed  of  hourly 
surface  wind  observations  (Lester  and  Phillips,  1987).  It 
follows  that  some  assumptions  must  be  made  about  the  nature 
of  the  structure  of  the  flow  aloft  in  order  to  carry  out  the 
proposed  experiment  for  the  Alberta  SDWs. 

Wind  profiles  have  been  specified  by  measuring  a 
'regionally'  representative  geostrophic  wind  at  700  mb  (near 
mountain  top  level)  and  then  linearly  interpolating  between 
that  wind  and  the  surface  wind.  The  requirements  for  the 
stability  profile  were  simplified  by  suppressing  the  criti¬ 
cal  streamline  feature  of  COMPLEX  and  ensuring  that  the  flow 
follows  the  shape  of  the  upper  boundary  which  was  predeter¬ 
mined  on  the  basis  of  mesoscale  flow  patterns  observed  else¬ 
where  during  SOW  conditions. 

The  slope  of  the  boundary  layer  top  for  COMPLEX  and  the 
final  flow-following  coordinate  system  was  specified  with  a 
simple  diagnostic  model  of  the  mountain  lee  wave  influence. 
Because  of  the  high  frequency  of  the  Chinook  arch  cloud  in 
southern  Alberta  (Chapter  2) ,  the  model  is  based  on  the 
assumption  that  the  SDW  cases  are  driven  by  a  hydrostatic 
lee  wave,  i.e.,  the  windstorms  are  primarily  the  result  of 
processes  described  by  Klemp  and  Lilly  (1975,  1978). 
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The  basic  characteristics  of  the  hydrostatic  wave 
(i.e.,  wavelength,  amplitude!  were  determined  from  past  sim¬ 
ulations  and  analyses  of  lee  wave  systems  associated  with 
downslope  windstorms  (e.g.,  Lilly  and  Zipser,  1972;  Finger- 
hut  and  Lester,  1973;  Klemp  and  Lilly,  1975,  1973;  Lester 
and  MacPherson,  1977).  Because  of  the  detailed  documenta¬ 
tion  and  analysis  of  the  11  January  1972  SOW  in  Boulder, 
Colorado,  and  because  that  case  apparently  captures  the 
major  mesoscale  wave  features  of  the  SDW,  it  was  subse¬ 
quently  chosen  as  the  primary  guide  for  configuring  the 
boundary  layer  top  of  COMPLEX.  The  lower  portion  of  the 
potential  temperature  field  for  the  11  January  1972  case, 
redrawn  from.  Lilly  and  Zipser  (1972)  is  presented  in  Fig. 

2a.  On  the  basis  of  the  spatial  variation  of  the  lowest 
stable  layer  in  that  analysis,  the  functions  illustrated  in 
Figure  2b  and  discussed  below  were  used  to  approximate  a 
two-dimensional,  ' hydros tacic  wave'  disturbance  of  the 
boundary  layer  top.  This  may  be  interpreted  as  a  condition 
in  which  the  model  flow  surfaces  correspond  to  isentropic 
surfaces,  and  where  small  scale  turbulent  mixing  plays  a 
negligible  role.  The  condition  is  introduced  via: 

Ttop(r)  “  htop  +  ztop  "  A  exp(-b[r  -  *tQpl2)  .  CO) 

A  =  ^top  +  ztop  ”  zbase  ~  ^base  '  Cl) 

where  Tcop(r)  is  the  height  (ASL)  of  the  boundary  layer  top 
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(a)  POTENTIAL  TEMPERATURE  ANALYSIS 

HciqrtC  ASL  (al 


(b'.i  TUFjGRAPHY  Ur  FLOW  SURFACES 

Meighc  ASL  (  r» ) 


Qljcjncr  East  (<•.169  km  per  grid  interval) 


Figure  2.  tj)  Potential  temperature  analysis  from  3oulder  windstorm  of 
11  January  1972  (after  Lilly  and  Zipser,  1972).  Potential 
temperature  isopleths  are  drawn  at  2.5  K  intervals.  The 
analysis  extends  from  the  Continental  Divide  (CD)  through 
Boulder  (Bou),  Jefferson  County  Airport  (JCA),  and  Denver 
(Den).  (b)  Schematic  cross  section  of  model  flow  surfaces 
illustrating  terrain  height  and  layer  thickness  parameters. 
The  upstream  tilt  of  the  hydrostatic  wave  trough  axis  is 
represented  by  a  thin  dashed  line. 
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at  r;  r  is  the  normal  distance  from  the  wave  trough  axis 
(positive  east  of  the  trough  axis);  hcop,  h^ase  are  the 
terrain  heights  ASL  at  the  top  and  base  of  the  slope;  Ztop, 
zbase  are  t*ie  boundary  layer  thicknesses  at  the  top  and  base 
of  the  slope;  b  is  a  constant  determining  the  half-width 
(wavelength)  of  the  wave  trough;  A  is  the  wave  amplitude; 
and  s>-3D  is  a  phase  distance  constant  included  so  that  an 
upstream  tilt  in  the  flow  surfaces  can  be  specified.  The 
Gaussian  exponential  function  is  chosen  because  of  its 
symmetry.  Furthermore,  the  symmetry  point  can  be  used  to 
represent  the  mountain  wave  trough,  and  the  half -width 
parameter  (bi  can  be  easily  adjusted  to  produce  disturbances 
of  different  wavelengths.  Finally,  wave  phase  can  be 
varied  as  a  function  of  height  by  adjusting  the  position  of 
the  symmetry  point  for  each  flow  surface. 

The  lowest  flow  surface  was  set  at  the  lowest  terrain 
elevation  within  the  model  domain.  Flow  surfaces  between 
the  boundary  layer  top  and  the  lowest  terrain  elevation  were 
specified  using: 


Ei 

±. 


T1(r)  =  Ei  -  A,  exp(-bCr  -  0^]“) 

htop  +  ztop  -  U(i-1)/(M-1)]  r  <  0 

fhtop  +  ztop  -  hminJ  Ci-1] / CN-1]  r  >  0 

A  r  <  0 

A  -  [E.:  (r  <  0)  -  Ei  (r  >  0)  ]  r  >  0 

■L  ± 


(12; 


(13) 


(14) 
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where  (  )  specifies  the  flow  surface  level;  Tj.(r)  is  the 
height  of  the  ith  flow  surface  ASL;  is  the  equilibrium 

height  of  the  ith  flow  surface;  Aa  is  the  amplitude  of  the 
ith  flow  surface;  hmj_n  is  the  elevation  of  the  lowest 
terrain  point  in  the  model  domain;  0 x  is  the  phase  distance 
constant  for  the  ith  flow  surface;  N  is  the  number  of  flow 
surface  levels  used  in  the  model  (Fig.  2). 

Values  of  boundary  layer  thickness  at  the  top  and  base 
cf  the  terrain  slope  ZZ0D,  2base  ^-ciriS-  (10),  (11),  (13)), 

were  alsc  derived  from  the  11  January  1372  Boulder  potential 
temperature  analysis  (Fig.  2)  as  a  first  approximation  of 
the  chinock  mountain  wave  flow.  Terrain  heights  at  the  top 
and  bottom  cf  the  slope,  ( hzop,  hbase)  were  derived  from  the 
terrain  data  set  for  southern  Alberta  in  order  to  scale  the 
wave  amplitude  to  the  model  terrain.  An  13-20°  per  km  up¬ 
stream  tilt  (phase  change  with  height)  was  specified  in  the 
model  flow  surfaces  based  on  a  simulation  of  the  11  January 
1972  Bculier  windstorm  by  the  Klemp  and  Lilly  (1973).  In 
order  to  position  the  mountain  wave  disturbance,  the  hori¬ 
zontal  axis  of  the  hydrostatic  wave  trough  was  aligned  par¬ 
allel  to  the  Continental  Divide  (337°  to  157°  true) ,  coinci¬ 
dent  with  the  strongest  hourly  averaged  surface  winds  (this 
is  in  agreement  with  observation  and  theory).  Wave  propaga¬ 
tion  was  assumed  to  occur  away  from  the  mountains  at  an 
angle  normal  to  the  Continental  Divide  with  the  wave  ampli¬ 
tude  remaining  constant.  Space-time  analyses  of  hourly 
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averaged  wind  speeds  along  the  Banff -Calgary  corridor  were 
then  used  to  approximate  the  speed  of  the  wave  trough. 

In  summary,  the  COMPLEX  wind  analysis  code  has  been 
modified  by  using  observed  ard  simulated  characteristics  of 
hydrostatic  mountain  waves  to  specify  the  combination  of 
mesoscale  forcing  of  the  surface  flow  features  and  the  mass- 
consistent  adjustments  intended  to  account  for  smaller  scale 
flow  features  produced  by  terrain, 
b .  Case  and  Data  Selection 

Two  chincok  SOW  cases,  were  selected  from  t.oe  Chinook 
Project  data  base:  3  November  1975  and  27  March  1936. 
Selection  criteria  were  (i)  the  occurrence  of  record  or  near 
record  surface  winds  over  southwestern  Alberta,  and  (ii) 
the  availability  of  a  relatively  large  data  set  (Lester  and 
Phillips,  1937;  .  The  defined  study  area  (Fig.  3;  is  located 
in  the  lee  of  the  Canadian  Rocky  Mountains  near  Calgary.  45- 
120  km  east  of  the  Continental  Divide.  Upper  air  data 
included  radiosondes  from  the  northwestern  U.S.  and  south¬ 
western  Canada,  available  pilot  reports,  and  conventional 
GOES  satellite  imagery.  Hourly  surface  wind  reports  were 
gathered  from  a  variety  of  sources  including  conventional 
first-order  surface  stations  operated  by  the  Atmospheric 
Environment  Service  (AES),  surface  weather  stations  operated 
by  the  University  of  Calgary  (UC) ,  automated  weather 
stations  constructed  as  part  of  the  1988  Calgary  Winter 
Olympics,  and  air  quali ty/meteorological  data  monitoring 
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’ir»  3.  Topographical  chart  o£  southwestern  Alberta  (after  Lester, 
1978a'/.  The  study  area  is  indicated  by  heavy  dark  lines 
while  the  transect  line  is  indicated  by  a  dashed  Line. 


With  one 


sites  operated  by  natural  gas  processing  plants, 
exception,  data  from  all  of  the  above  sources  are  subjected 
to  high  quality  control  and  considered  to  be  reliable.  Site 
exposures  and  some  quality  control  procedures  for  data  from 
gas  plants  have  resulted  in  less  reliable  data  for  some 
locations.  Temporal  and  spatial  consistency  were  used  as 
the  primary  tool  to  accept  or  reject  those  data. 

Data  selected  as  model  inputs  from  the  data  base 
described  above  are  hourly  surface  wind  observations,  analy¬ 
ses  of  upper  air  flow  every  12  hours  and  terrain  elevations. 
For  the  majority  of  the  surface  wind  stations,  only  eight 
wind  directions  (N,  NW,  W,  SW,  S,  SE,  E,  NE)  were  recorded. 
The  average  density  of  surface  stations  is  one  per  850  km* 
for  the  3  November  1975  case,  and  one  per  600  km2  for  the  27 
March  1986  case.  These  may  be  compared  to  the  average 
density  of  one  per  12,600  km2  for  conventional  weather 
stations  representing  a  15  fold  density  improvement  for  the 
3  Novemoer  1975  case,  and  a  21  fold  density  improvement  for 
the  27  March  1986  case.  Selected  stations  for  each  case  are 
described  below. 

c .  General  Meteorological  Conditions 

( 1 )  3  November  1975 

Synoptic  patterns  from  surface  and  700  mb  analyses 
(Fig.  4)  for  3  November  1975  at  1400  MST  and  1700  MST, 
respectively,  shew  an  intense  cross-mountain  gradient  in 
mean  sea  level  pressure,  relatively  warm  air  aloft  (indi- 
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Figure  4.  Synoptic  analyses  (a)  surface  analysis  at  1400  MST 
and  (b)  700  mb  analysis  at  1700  MST  for  3  November 
1975.  The  province  of  Alberta  and  the  study  area 
are  indicated  by  heavy  dark  lines.  Both  (a)  and 
(b)  incorporate  standard  plotting  models  and  conven¬ 
tional  labels.  In  (a),  solid  lines  are  isobars 
at  4  mb  intervals.  The  solid  line  with  triangular 
pips  is  a  cold  front.  In  (b),  solid  lines  are 
contours  of  geopotential  height  at  30  m  intervals 
while  dashed  lines  are  isotherms  at  5  C  intervals. 


32 


caced  by  ridging  of  isotherms  on  700  mb  analysis)  ,  and  a 
strong  westerly  upper  air  current  during  the  SDW  event. 
Previous  analyses  by  Lester  (1978a),  Lester  and  Bach  (1986), 
and  Lester  and  Phillips  (1987)  reported  that  the  strongest 
winds  occurred  with  the  approach  of  a  rapidly  moving  baro- 
clinic  disturbance  from  the  west.  The  air  mass  across 
southern  Alberta  was  unusually  warm  for  the  season,  as 
evidenced  by  a  high  tropopause  and  the  location  of  the  major 
jet  stream  axis  well  north  of  the  windstorm  area  (Lester  and 
Phillips  ,  1987)  . 

The  Chinook  arch  was  present  on  the  east  side  of  the 
Continental  Divide  over  the  windstorm  area  indicating  a 
hydrostatic  lee  wave  of  50-100  km  wavelength  (Lester  and 
Bach,  1936;  Lester  and  Phillips,  1987).  Also  present  during 
the  windstorm  were  shorter  lee  waves  (28  km)  embedded  within 
the  Chinook  arch.  The  shorter  waves  were  identified  with 
both  aircraft  data  and  satellite  observations  (Lester  and 
Bach,  138  6  1  . 

Record  surface  winds,  extensive  property  damage,  and 
forest  blowdown  occurred  as  the  windstorm  spread  eastward 
from  the  mountains  (Lester,  1978a).  Movement  of  the  wind¬ 
storm  at  the  longitude  of  Calgary  is  illustrated  in  a  space- 
time  section  of  hourly  averaged  wind  speed  (Fig.  5)  between 
Marmot  Creek  Basin  (MCB),  Kananaskis  Research  Centre  (ESC), 
Jumping  Pound  (JPd) ,  and  Calgary  International  Airport 
(YYC) .  The  space-time  sections  were  produced  using  spline 
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Time  (MST) 


»  I  i  1 

MC3  ESC  JPd  YYC 

Figure  5.  Space-time  section  of  mean  hourly  wind  speed  (m/s) 

for  3  November  1975  along  transect  line  between  Marmot 
Creek  basin  (MC3),  Kananaskis  Research  Centre  (ESC), 
Jumping  Pound  (JPd),  and  Calgary  International  Airport 
(YYC) Isotachs  are  drawn  at  5  m  s'  intervals,  except 
25  m  s  is  not  shown  for  clarity. 
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interpolation  and  applying  a  single  smoothing.  The  transect 
line  is  represented  by  the  dashed  line  in  Fig.  3.  Wind 
speeds  exceeding  15  m  s“^  appeared  at  MCB  at  0600  MST.  The 
windstorm  then  proceeded  eastward  reaching  its  maximum 
strength  with  winds  exceeding  25  m  s-^-  when  it  arrived  at 
JPd  at  1300  MST  (note:  for  clarity,  the  25  ms"^  isotach  is 
not  shown  in  Fig.  5).  The  windstorm  reached  midway  between 
JPd  and  YYC  at  1400  MST  producing  winds  exceeding  15  m  s-^, 
but  by  the  time  it  reached  YYC  at  1500  MST,  the  windstorm 
diminished  in  intensity  as  wind  speed  fell  below  15  m  s-~. 

(2)  27  March  1986 

Synoptic  patterns  from  surface  and  700  mb  analyses 
(Fig.  6;  for  27  March  1936  at  1400  MST  and  1700  MST  respec¬ 
tively  are  similar  to  those  from  3  November  1975.  Again,  an 
intense  cross -mountain  mean  sea  level  pressure  gradient, 
relatively  warm  air  aloft,  and  a  strong  westerly  upper  air 
current  were  present  during  the  SDW  event.  These  conditions 
were  conducive  to  the  production  of  mountain  waves,  and,  in 
fact,  Lester  and  Phillips  (1937)  documented  the  presence  of 
a  Chinook  arch  over  the  windstorm  area. 

The  space-time  section  analysis  of  hourly  averaged  wind 
speed  between  Fortress  Mountain  (FtM)  ,  ESC,  JPd,  and  YYC  is 
shown  in  Fig.  7.  Unlike  the  3  November  1975  case,  wind 
speeds  at  the  higher  elevations  of  FtM  and  ESC  were  less 
than  15  m  s-^,  and  stronger  winds  at  lower  elevations  were 
observed.  The  windstorm  produced  winds  exceeding  30  m  s~^ 
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JPd 


YYC 


Space-time  section  of  mean  hourLy  wind 
March  1986  aLong  transect  line  between 
(FtM),  Kananaskis  Research  Centre  (ESC) 
and  Calgary  International  Airport  (YYC) 
drawn  at  5  m/s  intervals. 


speed  (m/s)  for  27 
Fortress  Mountain 
,  Jumping  Pound  (JPd), 
Isotachs  are 
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when  it  arrived  at  JPd  at  1300  MST.  Winds  exceeding 
20  m  s'1  continued  as  the  windstorm  reached  midway  between 
JPd  and  YYC  at  1400  MST.  Winds  generally  diminished  as  the 
windstorm  reached  YYC  at  1600  MST,  although  wind  speed  still 
exceeded  15  m  s-^-. 

d .  Application  of  COMPLEX  to  SDW  Cases 

( 1 )  Computational  Grid,  Terrain,  and  Station  Locations 
In  the  present  study,  a  three-dimensional  grid 
(34x23x6)  with  a  horizontal  resolution  of  4.169km  was  speci¬ 
fied  for  the  Calgary  area  for  both  the  3  November  1975  case 
and  the  27  March  1936  case.  The  grid  interval  was  selected 
because  the  available  terrain  data  were  coded  in  minutes  of 
latitude/ longitude  and  COMPLEX  requires  that  terrain  data  be 
arranged  in  square  grid  elements. 

A  three-dimensional  perspective  of  the  study  area  ter¬ 
rain  is  presented  in  Fig.  8.  The  steepest  terrain  near  the 
Continental  Divide  was  not  included  in  the  computational 
domain  to  avoid  areas  with  a  high  probability  of  non-hydro- 
static  effects.  Important  terrain  features  in  Fig.  8  are 
the  Bow  River  Valley,  and  the  relatively  steep  secondary 
slopes  lying  west  of  Highway  2  near  Crossfield.  Maps  of  the 
model  grid  and  the  surface  wind  stations  for  each  case  are 
presented  in  Fig.  9.  The  eastern  portion  of  the  grid  is  not 
shown  because  no  surface  wind  data  were  available.  The 
names  and  grid  coordinates  of  each  surface  wind  station  are 
listed  in  Table  1. 
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Table  1.  Surface  wind  stations 
(a)  3  November  1975  case 

Station  Name  (Symbol)  Latitude (N)  Longitude (W) 


E.  Crossfield  (3) 

51.40 

114.07 

Carstairs  (14) 

51 . 57 

114.19 

Lone  Pine  Creek  (9) 

51.43 

113 . 37 

Harmattan  (7) 

51.67 

114.44 

Wildcat  Hills  (12) 

51.23 

114.63 

Jumping  Pound  NE  (4a) 

51.15 

114.55 

Jumping  Pound  MW  (4b) 

51.16 

114 . 59 

Jumping  Pound  SE  (4c) 

51.13 

114.55 

Univ.  of  Calgary  ( UC ) 

51.09 

114 . 17 

Calgary  Inti  Airport  (YYC) 

51 . 12 

114.02 

2"  March  1936  case 

Station  Name  (Symbol) 

La 

titude (N; 

Longitude (W) 

Calgary  Inti  A.  ('HO 

51.12 

114.02 

Crossfieid  #2  (2a, 

c '  ->2 

113.92 

Crossfieid  #3  (2b) 

51.13 

113.92 

Crossfield  #4  trailer 

(la) 

51 .40 

114.03 

E.  Crossfieid  (51! 

51.33 

114.01 

W.  Crossfieid  (52) 

51.33 

114.12 

Jumping  Pound  NE  (3c; 

51.15 

114.55 

Jumping  Pound  NW  (3d) 

51 .16 

114 . 59 

Jumping  Pound  SE  ( 3 e ) 

51.13 

114.55 

Lone  Pine  Creek  (10a) 

51.59 

113 .33 

Lone  Pine  Creek  S.  (5a 

) 

51.51 

113.36 

Springbank  A.  (33) 

51 . 10 

114 . 37 

Harmaotan  (7a) 

51.6" 

114 . 44 

Wildcat  Hills  (11a) 

51  22 

114.63 
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Observations  of  hourly  averaged  surface  winds  are  shown 
in  Figs.  10-13  for  the  two  cases.  Summaries  of  surface  and 
700  mb  wind  observations  are  presented  in  Table  2.  As  noted 
in  the  previous  section,  synoptic  patterns  for  the  two 
storms  are  similar,  however,  there  are  significant  differ¬ 
ences  in  the  wind  distributions.  Inspections  of  plocted  and 
tabulated  winds  in  Figs.  10-13  and  Table  2  shew  that  in  the 
3  November  197  5  case,  the  7GC  mb  wind  speeds  are  about  25% 
higher  than  in  the  2~  March  1936  case,  yet  the  observed 
surface  winds  are  higher  in  the  latter  case.  this  differ¬ 
ence  is  mere  apparent  when  only  surface  wind  speeds  from 
these  stations  common  to  both  cases  are  compared  (Table  2). 
for  example,  at  14CC  MOT,  the  root  mean  square  wind  speed 
for  those  stations  is  15.3  ms'"-  for  the  1975  case  and  19.2 
ms'*  for  the  1936  case.  In  terms  of  kinetic  energy  of  the 
observed  winds ,  the  latter  case  is  nearly  4C%  stronger. 

(  2  Model  Sensitivity  Tests 

An  area  of  great  uncertainty  in  the  application  of 
C3MP1EX  to  SOW  flow  is  the  relationship  of  the  lowest  model 
flow  surfaces  to  the  terrain.  For  situations  in  which  the 
flow  surfaces  intersect  terrain,  flow  proceeds  around  ter¬ 
rain  obstacles  rather  than  over  them.  One  problem  is  that 
the  degree  to  which  this  intersection  occurs  cannot  be 
determined  from  the  observational  data. 
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Table  2.  Observed  surface  and  700  mb  geostrophic  wind 
speeds  and  directions. 

(a)  3  November  1975 

Mind  direction (deg) /speed {m  s-^) 


700  mb  geostrophic  wind 

Station  Name  (Symbol) 

E .  Crossfield  (  3  > 

Carstairs  (14) 

Lone  Pine  Creek  (9) 
Harmatran  (7) 

Wildcat  Hills  (12) 

Jumping  Pound  NE  (4a) 
Jumping  Pound  NW  (4b) 
Jumping  Pound  SE  (4c) 

LTniv.  of  Calgary  (UC) 
Calgary  Inti  Airport  (Y YC) 


(bi  27  March  1936 


700  mb  geostrophic  wind 

Station  Name  (Symbol! 
Jumping  Pound  NE  (3c) 
Jumping  Pound  NW  (3d) 
Jumping  Pound  SE  (  'e) 

Lone  Pine  Creek  (10a) 
Crossfield  #2  (2a) 
Crossfield  #3  (2b) 
Crossfield  #4  trailer  (la) 
Lone  Pine  Creek  S.  (5a) 
Calgary  Inti  Airport  (YYC) 
E.  Crossfield  (51) 

W .  Crossfield  (52) 
Springbank  Airport  (33) 
Harmattan  (7a) 

Wildcat  Hills  (11a) 


1300  MST 

1400  MST 

1500  MST 

269/25.0 

263/25.0 

269/25.0 

270/18 . 3 

270/17.8 

270/17.3 

225/11.1 

225/17.4 

225/20 . 1 

130/11.1 

130/3.9 

130/3.0 

225/9 . 3 

225/12.0 

225/10.7 

225/17.8 

225/12.9 

225/11.6 

225/25.0 

225/17.3 

225/13 . 3 

130/23.2 

130/16 . 5 

130/14.3 

N/A 

225/25.0 

225/16 . 5 

270/12.0 

270/14 . 3 

270/13 . 4 

225/3 . 0 

253/12.9 

253/12,9 

i.nd  direct 

ion (deg) /speed (m  s_1 

1300  MST 

1400  MST 

1500  MST 

266/13.6 

267/13.7 

269/13.9 

225/33 . 6 

N/A 

N/A 

225/27.5 

N/A 

N/A 

225/34.7 

225/33.6 

225/22.5 

270/13.0 

270/14 . 4 

270/14 . 4 

270/15 . 2 

225/16.1 

270/13. 3 

270/13.3 

270/15 . 2 

270/13 . 0 

130/13.3 

130/16 . 3 

130/15.0 

225/16.6 

225/16.6 

225/13 . 0 

270/11.3 

240/11.9 

250/16 . 3 

247/15.2 

247/17 . 5 

247/19.7 

247/13 . 3 

247/13 . 3 

247/20.0 

260/20 . 5 

250/20.5 

240/15.5 

270/6.9 

270/6.9 

270/6.9 

N/A 

N/A 

N/A 
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In  the  formulation  of  the  flow  surface  configuration 
{Eqns.  10-14),  the  thicknesses  between  the  model  flow  sur¬ 
faces  and  the  terrain  surface  are  primarily  controlled  by 
the  average  boundary  layer  thickness,  Z,  and  the  exponential 
constant,  b  (Eqn.  (12)),  so  model  sensitivity  experiments 
were  performed  to  test  both  the  effect  of  varying  layer 
thicknesses  (Sensitivity  tests  #1,  #2,  and  #3)  and  varying 
wavelength  (Sensitivity  tests  #4,  and  #5)  for  the  time  when 
the  windstorm  was  at  JPd.  The  27  March  1936  case  was 
selected  for  the  tests  because  it  contains  the  most  surface 
wind  observations. 

An  additional  problem  arose  in  the  formulation  of  the 
flow  surface  configuration.  Because  the  wave  trough  is 
forced  to  progress  eastward  without  changing  its  shape,  the 
flow  surfaces  rise  as  much  as  600  m  (over  50%  of  the  average 
boundary  layer  thickness)  near  the  base  of  the  terrain 
slope.  This  behavior  does  not  agree  with  observations 
(e.g.,  Lilly  and  Zipser,  1972).  In  order  to  present  a  more 
realistic  surface  wind  field  upstream  of  the  wave  trough. 
Ludwig  (1937)  suggested  that  the  wave  shape  be  modified  for 
cases  in  which  the  wave  trough  is  located  through  Btw  and 
YTC.  As  a  first  approximation,  the  value  of  b  upstream  of 
the  wave  trough  was  adjusted  so  that  no  change  in  flow  sur¬ 
face  height  occurred  at  the  southeast  corner  (anchor  point) 
of  the  grid  as  the  wave  trough  was  progressed  eastward 
toward  Btw  and  YYC . 
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Table  3  lists  function  parameters  held  constant  for  all 
runs,  and  their  values.  Values  of  parameters  that  were 
changed  for  sensitivity  tests  #l-#5  are  listed  in  Table  4. 

Final  analyses  of  the  surface  wind  fields  with  COMPLEX 
were  accomplished  for  each  case  at  three  separate  times: 

(i)  when  the  windstorm  was  located  at  JPd  (Fig.  3),  (ii) 
when  the  storm  reached  a  point  midway  between  JPd  and  YYC 
(3tw)  ,  and  (iii)  when  it  reached  Y“YC.  In  subsequent  sec¬ 
tions,  these  analyses  will  be  referred  to  as  Runs  #1,  #3, 
and  #3,  respectively.  As  described  earlier,  the  movement  of 
the  windstorm  was  simulated  by  placing  the  wave  trough  axis 
parallel  to  the  Continental  Divide  and  through  each  location 
noted  above.  The  space-time  analyses  of  hourly  averaged 
wind  speed  shown  in  Figs.  5  and  7  were  used  to  approximate 
the  time  at  which  the  windstorm  arrived  at  these  locations. 
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Table  3.  Mountain  wave  flow  parameters  held  constant  for 
all  model  runs.  Values  were  derived  from  Lilly 
and  Zipser  (1972)  potential  temperature  analysis 
of  11  January  1972  Boulder  windstorm. 


Parameter 

Value 

""top 

1270  m 

zbase 

762  m 

^top 

2134  m 

1ad  ase 

1219  m 

Table  4.  Mountain  wave  parane 
sensitivity  tests 

ters  specified 

in  mode 

Sensitivity  Effective 

Test  No.  Wavelength  (km) 

Paramete 
b  ( km  “  * ) 

r- 

Z  (n) 

1 .  75  . 

.0013 

1270  . 

2 .  75  . 

.0013 

1000  . 

3  .  75  . 

.0013 

1100  . 

4  .  50  . 

.0044 

1100  . 

2  .  inn 

.0011 

1100  . 
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Chapter  4 


RESULTS 


a .  Sensitivity  Tests 

{ 1 }  Effect  of  Boundary  Layer  Thickness 
In  tests  #1,  #2,  and  #3,  the  value  of  b  in  Eqn.  16  was 
set  to  produce  an  effective  wavelength  of  75  km  and  the 
boundary  layer  thickness  was  systematically  varied.  East- 
west  cross  sections  of  the  model  flow  surfaces  and  the 
terrain  surface  through  JPd  and  YYC  are  presented  in  Fig.  14 
for  tests  #1  and  #2.  The  respective  surface  wind  vectors 
and  isotach  analyses  are  shown  in  Figs.  15  and  16.  In  each 
analysis,  the  wave  trough  axis  is  represented  by  a  heavy 
dashed  line.  For  clarity,  surface  wind  vectors  are  dis¬ 
played  at  8.333  km  intervals  (odd  numbered  rows  and 

_  * 

columns)  .  Isocachs  are  plotted  at  a  5  m  s  *  intervals  and 
areas  of  wind  speeds  exceeding  15  a  s”^-  are  cross-hatched . 

In  test  #1  Z  was  set  equal  to  (1270  m)  to  deter¬ 

mine  the  effect  of  allowing  the  flow  surfaces  to  intersect 
the  terrain.  Wind  vectors  in  Figure  15  indicate  flow  is 
increased  near  the  wave  trough  axis,  particularly  through 
the  Bow  River  Valley  near  Wildcat  Hills  (station  11a)  where 
westerly  wind  speeds  exceed  25  m  s~^.  The  largest  gradients 
in  wind  direction  are  noted  in  the  Bow  River  Valley  and 
along  the  northern,  western,  and  southern  boundaries.  A 
speed  minimum  {<_  10  m  s-^)  is  found  over  the  hilltops  near 
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TOPOGRAPHY  OF  FLOW  SURFACES 

Height  ASL  (m) 


TOPOGRAPHY  OF  FLOW  SURFACES 


Figure  14.  East-west  vertical  cross  section  through  Jumping  Pound 
(JPd)  and  Calgary  International  Airport  (YYC)  of  model 
flow  surfaces  for  (a)  Sensitivity  test  #1  and  (b) 
Sensitivity  test  §2.  Lines  from  bottom  to  top  are 
the  terrain  surface  and  the  five  upper  r,odel  flow 
surfaces.  Note:  Sensitivity  test  #3  produced  optimum 
results.  The  corresponding  cross  section  is  presented 
just  prior  to  the  case  studies  (Figure  22). 
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Same  as  Fig.  If.  except  for  Sensicivity 


51.40N, 


occur 


114. 30W  while  wind  speeds  exceeding  10  m  s 
between  that  location  and  Highway  2  to  the  east. 

In  general,  test  #1  shows  several  realistic  results  due 
to  both  mesoscale  forcing  and  topography:  (i)  increased 
winds  are  near  the  wave  trough,  ( i i )  WSW  wind  directions 
with  small  scale  variations  due  to  terrain  features,  and 
(iii)  numerous  wind  speed  maxima /minima  which  reflect 
terrain  details.  However,  medal  wind  speeds  substantially 
underestimate  observed  wind  speed  values  (Fig.  10,  Table  2) 
at  J?d  (stations  3c,  3d,  and  3e),  Spri.ngbank  Airport 

(station  33) .  and  V.  Crossfield  (station  52)  bv  30-75i 
(4-21  xs“^!  . 

In  test  #2,  Z  was  reduced  to  1300  m  to  reduce  the 
degree  of  terrain  intersection  of  the  flow  surfaces  at  the 
base  of  the  terrain  slope.  Patterns  of  model  winds  relative 
to  wave  trough  location  are  similar  to  these  produced  in 
test  #1  (Fig.  15)  except  that  fewer  winds  exceed  20  m  s'", 
and  there  is  a  larger  area  of  winds  exceeding  15ms  *  along 
the  wave  trough  axis  near  J?d.  Wind  direction  varies  little 
through  the  grid  domain  except  at  the  boundaries  and  over 
the  hills  just  west  of  Crossfield  # 2  (Station  2a;  .  Model 
wind  speeds  near  JPd  (17  m  s_ and  Springbank  (12  m  s  ^ ) 
are  closer  to  the  observed  values  then  in  test  #1.  Greater 
wind  speeds  also  occur  at  the  base  of  the  mountains  near  the 
southwest  grid  corner.  Problems  still  exist  near  W.  Cross- 
fieid  where  the  model  wind  speed  underestimates  the  observed 
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wind  speed  even  more  than  in  test  #1,  and  the  maximum  model 
wind  speed  severely  underestimates  observed  maximum  wind 
speed  (20  vs  34  m  s-^). 

In  test  #3,  Z  was  set  to  1100  m,  i.e.,  an  intermediate 
value  between  tests  #1  and  #2.  This  test  produced  the  best 
results.  The  west-east  distribution  of  model  wind  speed 
through  JPd,  Springbank  Airport  (Spr),  and  YYC  are  shown  in 
Fig.  I”  for  tests  #1,  #2,  and  #3.  Variation  in  wind  speed 

response  to  terrain  from  tast-to-test  is  apparent,  as  from 
Spr  to  YYC  wind  speed  increases  and  decreases  dc  not  occur 
at  tr.e  same  terrain  points.  Near  JPd,  ail  three  tests  pro¬ 
duce  wind  speed  maxima,  but  only  in  test  #3  is  the  maximum 
model  wind  speed  (30  r.s"*j  close  to  the  observed  value 
(34  m  s'1)  . 

( 2 '  Effect  of  Lee  Wave  Length 

In  tests  #4  and  #5.  the  value  of  2  was  fixed  at  1100  m 
while  the  value  of  b  was  varied  to  produce  effective  wave¬ 
lengths  cf  50  km  and  ICC  km.  respectively.  East-west  cross 
sections  of  the  model  flow  surfaces  through  JPd  and  YYC  are 
presented  for  both  tests  in  Fig.  13.  The  respective  surfac 
wind  vectors  and  isctach  analyses  are  presented  in  Figs.  19 
and  20  . 

In  test  #4  (Fig.  19),  the  relatively  short  wavelength 

produces  a  narrow  region  cf  increased  surface  wind  flow. 

Wind  speed  maxima  ( >_  20  m  s-^)  occur  in  line  with  the  wave 

—  1 

trough  axis.  Maximum  model  wind  speeds  (27  ms*)  occurs 
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near  51.25N,  114. 70W.  Downstream  of  the  wave  trough,  wind 
speed  and  direction  show  small  gradients.  Only  slight  vari¬ 
ations  in  wind  direction  occur  near  JPd,  while  greater  wind 
shifts  occur  along  the  lateral  boundaries. 

Several  poor  results  were  noted  in  test  #4.  These 
included:  (i)  maximum  model  wind  speed  (27  m  s-^)  are  less 
than  80%  of  the  observed  value  (34  m  s-^)  ;  (ii)  model  wind 
speeds  underestimate  observed  values  at  JPd,  Springbank,  and 
W.  Cross f i  i Id ;  and  (iii)  little  model  response  to  terrain 
features.  Only  near  the  wave  trough  axis  do  the  model  wind 
vectors  appear  to  respond  realistically  to  terrain  varia¬ 
tions  . 

The  longer  wavelength  in  test  #5  (Pig.  20),  produces 
increased  model  wind  speeds  downstream  of  the  wave  trough 
axis.  Patterns  of  wind  speed  maxima/minima  similar  to  Pig. 
15  are  produced,  except  that  high  winds  occur  south  of  JPd. 
Maximum  model  wind  speeds  (33  m  s-1)  are  found  at  51.05N, 
114. SOW.  In  addition,  more  realistic  variations  in  wind 
direction  occur  in  the  Bow  River  Valley  near  Wildcat  Hills 
and  down  the  steep  slope  from  the  southwest  corner  of  the 
grid . 

Generally,  test  #5  shows  that  an  increase  in  the  hori¬ 
zontal  scale  of  mesoscale  forcing  causes  an  increase  of  the 
horizontal  extent  of  strong  surface  winds.  In  addition,  the 
m  iximum  model  wind  speed  is  comparable  to  the  observed  value 
(33  vs  34  m  s-i-}  and  model  winds  are  more  responsive  to 


61 


terrain  features.  However,  some  of  the  problems  in  other 
tests  also  occur  in  test  #5.  Observed  speeds  near  JPd, 
Springbank,  and  W.  Crossfield  are  significantly  underesti¬ 
mated  . 

In  an  effort  to  position  the  wind  speed  maximum  closer 
to  JPd,  the  value  of  b  was  adjusted  to  produce  an  effective 
wavelength  of  75  km  with  Z  set  at  1100  m  (Test  #3)  .  The 
west-east  distribution  or  model  winds  through  JPd,  Spring- 
bank  (Spr) ,  and  YYC  are  presented  for  tests  #3,  #4,  and  #5 
in  F^g.  21.  All  three  tests  show  increases  in  wind  speed 
near  JPd.  but  only  test  42  produces  a  model  wind  speed  com¬ 
parable  to  the  observed  value  at  that  location.  Thus,  the 
combined  results  of  all  tests  (Figs.  17,  21)  shew  clearly 
that  test  #2  produces  the  best  results. 

In  the  two  case  studies  which  are  presented  in  the 
following  section,  flow  surface  parameters  (b, Z)  defined  by 
sensitivity  test  42  are  used  to  represent  the  windstorm  at 
JPd  (Run  42,  .  As  the  storm  "moves  eastward,"  i.e.,  as  the 
trough  is  subsequently  positioned  at  Btw  (Run  #2)  and  YYC 
(Run  42)  ,  flow  surface  shapes  were  modified  so  that  no 
change  in  flow  surface  height  occurred  at  the  southwest  grid 
corner.  The  resulting  function  parameters  for  each  case  run 
are  specified  in  Table  5  while  east- west  cross  sections  of 
the  model  flow  surfaces  through  JPd,  Btw, 
presented  in  Figs.  22  and  23. 

5Z 


and  YYC  are 


Table  5.  Mountain  wave  flow  parameters  specified  for 

3  November  1975  and  27  March  1986  case  studies 
(see  Figs.  22  and  23a, b). 


Case 
Run  No . 

Wave 

Position 

b  [upstream] 
<10~3  km-2) 

Parameter 

b [downstream] 
do-3  km"2) 

Z 

(m> 

Fiqure 

1 . 

JPd 

1.800 

1.8 

1.1 

22 

Btw 

0.767 

1 . 8 

1.1 

23a 

3  . 

YYC 

0 . 420 

1 . 3 

1 . 1 

23b 
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TOPOGRAPHY  OF  FLOW  SURFACES 

Heighc  ASL  (m) 


JPd  YVC 

Distance  East  (^.’69  km  per  grid  interval) 


East-west  vertical  cross  section  through  Jumping  Found  (JPa! 
and  Calgary  Ince rnac  iona 1  Airport  CffC)  of  model  flow  surf¬ 
aces  fjr  Case  Hun  #1  (Sensitivity  test  #3).  Lines  from 
bottom  to  top  are  the  terrain  surface  and  the  five  upDer 
model  flow  surfaces. 
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b. 


Case  Studies 


( 1 )  3  November  1975 

Model  anemometer  height  wind  vectors  and  isotach  analy¬ 
ses  are  presented  for  case  runs  #1,  #2,  and  #3  in  Figs.  24, 
25,  and  26,  respectively.  For  ease  in  comparison,  these 
figures  are  grouped  at  the  end  of  each  case  discussion.  In 
this  and  in  the  27  March  1986  case  study,  the  position  of 
the  wave  trough  is  indicated  by  a  heavy  dashed  line  and 
areas  of  wind  speeds  exceeding  15  m  s-*  are  indicated  by 
cross-hatching. 

Cn  3  November  1975  the  wave  trough  moved  across  the 
western  portion  of  the  domain  to  arrive  at  JPd  at  1300  MST 
(Fog.  24  ,  Btw  at  1400  MST  (Fig.  25/,  and  finally  YYC  at 
1500  MST  (Fig.  26)  .  As  ncted  in  the  discussion  of  Fig.  5. 
The  times  are  approximate;  Lester  (1973a)  provides  some 
evidence  which  suggests  that  the  windstorm  moved  more 
rapidly  as  it  moved  past  Calgary. 

As  expected,  the  distribution  of  wind  speeds  at  each 
time  is  roughly  symmetric  with  respect  to  the  position  of 
the  wave  trough.  As  the  wave  trough  is  moved  from  JPd  to 
YYC,  model  wind  speeds  are  stronger  near  the  trough.  The 
highest  wind  speeds  and  wind  speed  gradients  appear  upstream 
of  the  trough  axis.  Downstream  ,  away  from  the  wave  trough, 
wind  speeds  decrease  gradually.  Wind  directions  are  mostly 
SW-WSW  with  greatest  variations  along  the  wave  trough  near 


67 


relatively  steep  terrain  features  and  along  the  lateral 


boundaries . 

Wind  speed  maxima  and  minima  are  very  localized  and 
occur  near  both  local  maxima  and  minima  in  terrain  eleva¬ 
tions.  This  variability  is  well-illustrated  in  Figs.  25  and 
26.  In  Fig.  25,  a  wind  speed  maximum  (>_  20  m  s"^)  occurs 
directly  ever  a  terrain  peak  {>_  13G0  mi  at  51.05N,  114. 25W 
with  lighter  winds  in  the  lower  elevations  surrounding  the 
peak.  Near  a  higher  terrain  peak  (>_  1400  m)  at  51.22N, 
114.25W,  wind  speed  is  hi  jr.er  around  the  peak  than  it  is 
directly  ever  it.  Similarly,  in  Fig.  26,  a  wind  speed  maxi¬ 
mum  (  25  m  s~~<  octurs  at  a  relatively  lower  terrain  eleva¬ 
tion  between  two  hill  tops  i_  1400  mi  at  51.22N,  114. 30W 

with  relatively  lighter  winds  directly  ever  the  terrain 
pea.es.  This  may  be  an  artifact  that  occurs  near  intersec¬ 
tions  or  a  flow  surface  with  the  terrain.  Adjustment  toward 
r.on-divergence  may  lead  to  unrealistically  high  speeds  in 
terrain  gaps. 

The  hignest  model  wind  speeds  occur  with  the  wave 
positioned  t .trough  J?d  (Fig.  24  The  maximum  model  wind 
speed  ( 2d  n  s-*  occurs  at  52.3~N.  114. 50W  and  20  m  s'" 
winds  appear  at  51 . 1  IN,  114. 50W  near  J?i  vacations  4a,  4b, 
ar.d  4c;  . 

As  the  wave  trough  moved  eastward  to  Btw  (Fig.  25  . 

areas  at  winds  exceeding  20  m  s decrease.  Maximum  model 


rind  speed  (25  ms1)  with  the  trough  at  BTW  occurs  at 


51.17N,  113. 31W. 

When  the  wave  trough  is  moved  to  YYC  (Fig.  26),  the 

maximum  wind  speed  (26  m  s-1)  appears  at  51.31N,  114. 30W, 

within  a  larger  area  of  winds  exceeding  15  m  s_i.  Areas  of 

_  -1 

winds  exceeding  20  m  s  x  have  diminished  in  size  and  number. 

A  comparison  of  case  runs  #1,  #2,  and  #3  (Figs.  24,  25, 
and  25)  with  the  respective  observed  surface  winds  (Figs. 

10,  11a,  and  lib!  shows  that  the  model  successfully  reprc- 


high  winds  observed  at  the  base  of  the  major  terrain 


slice  for  the  case  in  which  the  wave  trough  is  located 
through  J?d.  Additionally,  model  winds  near  Harmattan 
•station  7 )  and  Carscairs  (station  14)  respond  realisti¬ 
cally,  i.e.,  they  increase  as  the  wave  trough  is  moved  east¬ 
ward.  seed  agreement  between  model  and  observed  wind  speeds 
is  obtained  near  YYC  for  ail  three  case  runs.  Underesti¬ 
mates  of  observed  wind  speeds  by  the  model  occur  at  J?d  for 
cases  in  which  one  wave  trough  is  located  to  the  east,  and 
at  Crossfield  (station  3 ;  for  all  three  runs.  In  each  of 
tnese  situations,  the  stations  are  toe  far  away  from  the 
wave  trough  to  be  affected  by  the  model  wave  influence. 

This  indicates  that  tne  actual  flow  structure  may  feature  a 
mountain  wave  mode  wnose  influence  (i.e.,  mesoscale  forcing 


of  boundary  layer  thickness)  is  horizontally  mere  extensive 


(2)  27  March  1986 


Model  anemometer  height  wind  vectors  and  isotach  anaiy 
ses  are  presented  for  case  runs  #1,  #2,  and  #3  in  Figs.  27, 
23,  and  29,  respectively.  The  wave  trough  was  moved  across 
the  western  portion  of  the  domain  to  arrive  at  JPd  at  1300 
MST  (Fig.  27 } ,  Btw  at  1400  MST  (Fig.  23),  and  finally  YYC  at 
1500  MST  (Fig.  29 )  . 


There  are  several  similarities  between  this  case  and 
the  3  November  1575  case.  The  surface  wind  speed  distribu¬ 
tion  at  each  time  is  symmetric  with  respect  to  the  position 
cf  tne  wave  trough.  As  the  wave  trough  moves  eastward  from 
JFd  to  YYC,  model  wind  speeds  increase  near  the  trough  and 
decrease  elsewhere.  The  highest  wind  speeds  and  wind  speed 
gradients  appear  upstream  of  the  trough  axis.  Also,  flew  is 
most ly  SW-W2W  with  the  greatest  variations  in  wind  direction 
along  the  wave  trough  near  relatively  steep  terrain  features 
and  along  the  lateral  boundaries.  Finally,  wind  speed 
maxima  and  minima  are  very  localized,  and  may  occur  near 
either  local  maxima  or  minima  in  terrain  elevation. 

A  di iterance  between  this  case  and  the  3  November  case 
appears  when  the  wave  is  located  through  both  JPd  and  Btw 
(Figs.  2~  and  23  .  Secondary  wind  maxima  (>  10  m  s 

appear  downstream  cf  the  wave  trough  axis.  Model  wind 
speeds  decrease  below  10  m  s-*  between  51.15N,  114. 30W  and 
5L.55N.  11 4. BOW  and  increase  above  10  m  s-^  immediately 

downstream.  The  boundaries  of  the  secondary  maxima  appear 
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to  follow  the  1200  m  and  1300  m  terrain  contours. 

As  in  the  3  November  case,  the  highest  model  wind 

speeds  are  associated  with  the  wave  position  through  JPd 

(Fig.  27) .  Maximum  model  wind  speed  (30  m  s'^)  occurs  at 

51. 10N,  114. SOW  near  JPd  (stations  3c,  3d,  and  3e). 

As  the  wave  trough  is  moved  eastward  to  Btw  (Fig.  2 S )  , 

_  1 

areas  of  winds  exceeding  20  ms*  decrease  in  size  and 


rr.be  r  . 

Maximum 

model  wind  speed  occurs 

at  51. 

.  17N, 

113 . 31 W 

eases  to 

24  m  s'*. 

As 

the  wave 

trough  reaches  YYC  (Fig. 

2  9  ) 

the 

maximum 

model  wind  speed  (25  m  s~*i  is  produced  at  51.31N,  114. 30W. 
The  area  of  winds  exceeding  15  m  s'*  is  larger  than  the  pre¬ 
vious  case  run,  but  areas  of  winds  exceeding  20  m  s'*  have 
further  decreased  in  size  and  number. 

A  comparison  of  case  runs  4 1,  #2,  and  #3  (Figs.  27,  28, 
and  23;  wic.n  the  respective  observed  surface  winds  (Figs. 

12,  13a,  and  13b:  shows  that  COMPLEX  again  accurately  repro¬ 

duces  the  hi gh  winds  observed  at  the  base  of  the  major 
terrain  slope  for  the  case  in  which  the  wave  trough  is 
located  through  JPd.  Although  winds  near  Harmattar.  (station 
7a)  are  also  accurately  reproduced,  underestimates  of 
observed  wind  speed  occurs  at  virtually  all  other  stations. 
In  particular,  model  wind  speeds  are  lower  than  observed 
winds  at  JPd  for  cases  when  the  wave  trough  is  located  to 
the  east,  i.e.,  at  Btw  and  T!C .  at  E.  Crossfield  (station 
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51),  Crossfield  #2  (station  2a),  and  Lone  Pine  Creek  S. 
(station  5a)  for  all  three  case  runs. 

(  3  )  Case  Comparison  and  Discussion 

A  comparison  of  model  output  winds  from  the  3  November 
1975  case  (Figs.  24-26)  with  those  from  the  27  March  1936 
case  (Figs.  27-29)  shews  that  greatest  differences  between 
the  two  cases  appear  downstream  of  the  wave  trough  when  it 
is  located  through  JPd  (Figs.  24  and  27)  and  Btw  (Figs.  25 
ar.d  23:  .  Downstream  of  the  wave  trough  in  the  27  March 
case,  winds  decrease  sharply  to  a  minima  of  <  10  m  s~~, 
while  in  the  3  November  1975  case,  the  decrease  in  wind 
speed  is  more  gradual-  Since  identical  flow  surface  config¬ 
urations  are  used  for  both  cases,  differences  between  each 
case  can  be  related  to  differences  in  the  input  surface  and 
upper  air  wind  data  (Figs.  10-13,  Table  2) .  For  the  27 
March  1936  case,  the  wind  speed  applied  at  the  boundary 
layer  top  is  6  m  s~~  less  than  that  of  the  3  November  1975 
case.  Additionally,  surface  wind  data  are  not  available 
from  Carstairs  (station  14)  and  surface  wind  speed  at 
Harmattan  (station  7a )  was  3-5  m  s  ~  ^  below  values  from  the  3 
November  1975  case.  Collectively,  in  the  27  March  1936 
case,  lighter  winds  at  the  boundary  layer  top,  and  the 
weighting  of  the  initial  wind  analysis  by  lighter  winds  at 
Harmattan  appear  to  be  the  cause  of  decreased  model  wind 
speeds  downstream  of  the  wave  trough. 
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Similarities  between  the  two  cases  appear  in  all  analy¬ 
ses  along  the  wave  trough  axis  and  when  the  wave  trough  is 
located  through  YYC.  For  both  cases,  wind  speed  maxima 
occur  at  nearly  the  same  locations  just  upstream  of  the  wave 
trough.  Locating  the  wave  trough  through  YYC  produces  only 
slightly  different  wind  speed  and  direction  pa  .terns  for 
both  cases.  Along  Highway  2  at  51.30N,  model  wind  speed  is 
higher  for  the  27  March  1936  case  while  at  51.60N,  114. 25W, 
model  wind  speed  is  higher  for  the  3  November  1975  case. 

This  difference  appears  to  be  related  to  differences  in  the 
input  surface  wind  data.  In  the  27  March  1936  case,  model 
winds  respond  to  inputs  from  Crossfield  #2  (station  2a)  and 
Crcssfield  #3  (station  2b.  which  are  not  available  for  the  3 
November  1975  case.  Similarly,  in  the  3  November  1975  case, 
model  winds  respond  to  inputs  from  Carstairs  (station  14) 
which  are  not  available  for  the  27  March  1936  case. 

In  general,  both  case  analyses  showed  increased  surface 
wind  flow  due  to  mesoscale  forcing  of  the  boundary  layer  by 
a  'hydrostatic  wave'  and  from  smaller  scale  topographic 
effects.  The  horizontal  extent  of  the  increased  model  sur¬ 
face  wind  flow  was  shown  to  vary  with  the  effective  wave¬ 
length  of  the  hydrostatic  wave.  Maximum  model  surface  wind 
speeds  occurred  at  the  bas*.  of  the  principal  terrain  slope 
and  were  within  6-12%  of  -the  maximum  observed  surface  wind 
speed . 
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Both  analyses  also  displayed  two  general  problems:  (i) 

there  were  unrealistic  changes  in  wind  direction  due  to  the 
lateral  boundaries  and  (ii)  there  was  considerable  under¬ 
estimation  of  observed  wind  speeds  upwind  of  the  wave  trough 
as  the  wave  progressed  from  JPd  to  YYC .  To  illustrate  the 
wind  speed  underestimation  problem,  two  west-east  distribu¬ 
tions  of  model  wind  speed  through  JPd,  Spr,  and  YYC  are 
presented  in  Fig.  30.  These  diagrams  are  constructed  from 
the  2"  March  1936  case  when  the  hydrostatic  wave  trough  axis 
was  located  through  Btw  (Fig.  30a)  and  YYC  (Fig.  30b) .  Both 
modified  (Fig.  30)  and  unmodified  (Chapter  3)  flow  surface 
shapes  were  examined,  however  model  wind  speeds  were  still 
toe  low,  i.e.,  the  flew  surface  shape  modification  used  in 
this  study  only  slightly  improves  the  model  wind  speed  up¬ 
wind  of  the  wave  trough  position. 

In  reference  to  the  latter  problem,  Durran  and  Klemp 
(193";  suggest  that  the  amplification  of  downslope  winds  is 
fundamentally  similar  to  the  amplification  produced  during 
the  transition  from  subcritical  to  supercritical  flow  in 
shallow-water  hydraulic  theory.  Results  from  their  non¬ 
hydrostatic,  two-dimensional  model  indicate  mesoscale  forc¬ 
ing  to  be  horizontally  more  extensive  than  that  produced  by 
a  hydrostatic  wave,  as  used  in  the  present  case.  It  is 
suggested  that  some  sort  of  similar,  'hydraulic  jump' 
forcing  would  eliminate  the  problems  with  underestimated 
wind  speeds  upstream  cf  Che  mesoscale  wave  trough  as  dis- 


80 


JPd  Spr  YYC  E 

Distance  East  (u.163  km  per  grid  incerval) 


East-west  v.irucion  of  model  anemometer  height  wind  speed 
through  Jumping  Pound  (JPd),  Spnngbank  (Spr),  and  Calgary 
internal  ion,'.  1  Airport  (YYC)  for  modified  (dashed  line)  and 
unmodified  (solid  line  flow  surface  configurations  upstream 
of  wave  troughs  located  at  (a)  Btw  and  (b)  YYC.  Observed 
wind  speeds  are  denoted  by  an  asterisk  (See  text  for 
discussion)  . 
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cussed  above.  Future  experiments  with  COMPLEX  may  benefit 
from  a  modification  of  the  boundary  layer  top  to  reflect  the 
hydraulic  jump  configuration. 


3: 


Chapter  5 


SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 
FOR  FUTURE  RESEARCH 


A  mass-c 
red: fied  to  a 
stem  cases, 
its  r '  fixation 
order  to  situ 
tain  waves. 


t  w  c  —  c  i  m.  e  n  s  i  o  n 
previously  or 
patterns .  Th 


onsistent,  diagnostic  wind  model  (COMPLEX)  was 
naiyze  flow  during  two  Alberta  Chinook  wind- 
The  principal  change  to  COMPLEX  involved  the 
of  its  flow-following  coordinate  system  in 
late  forcing  of  the  surface  wind  flow  by  mcun- 
Specifically,  the  forcing  was  provided  by  a 
al  Gaussian  exponential  function  derived  from 
served  and  simulated  hydrostatic  mountain  wave 
e  principal  function  variables  were  the  average 
r  thickness,  and  effective  dis romance  wave- 
mum  values  of  the  parameters  were  determined 
of  sensitivity  tests. 


the  final  analyses  of  the  available  surface  wind  data 
were  ac r crop lisned  by  orienting  the  trough  line  of  the  two- 
dimensional  hydrostatic  wave  parallel  to  the  mountain  crest 
at  t.ne  point  of  the  maximum  observed  surface  winds.  The 
analyses  were  then  performed  for  each  case  at  hourly  inter¬ 
vals  between  the  foothills  and  Calgary.  Input  data  con¬ 
sisted  of  hourly  averaged  surface  winds  and  700  mb 


geostropnio  winds  . 

In  general,  objective  analysis  with  COMPLEX  produced 
realistic  wind  distributions  as  the  Chinook  windstorms 
progressed  towards  Calgary.  Maximum  winds  during  the  wind- 
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storms  ware  observed  at  the  foot  of  the  mourtains  and 


decreased  as  the  windstorms  moved  eastward.  Although  as 
would  be  expected  a  COMPLEX  su-  face  isotach  pattern  arose 
due  to  terrain  variations,  the  overall  influence  of  the  two- 
dimensional  "hydrostatic"  wave  was  clearly  evident  in  tne 
roughly  symmetrical  isotach  pattern  about  the  maximum  winds 
near  meso-scale  wave  trough. 

A  problem  cf  unrealistic  inflow  and  outflow  at  the 
lateral  boundaries  was  also  noted  where  the  lower  flow  sur¬ 
face  intersected  the  terrain.  This  problem  can  be  easily 
minimised  by  keeping  the  model  domain  significantly  larger 
t h a n  the  area  of  interest. 


Tne  observed  wind  analysed  with  COMPLEX  were  hourly 
average  winds  in  the  main,  i.e.,  no  consideration  was  giver, 
to  gusts  which  may  tc  exceed  twice  tha  ncurly  mean  speed 
especially  ir.  rugged  terrain.  In  fact,  the  inability  of 
COMPLEX  tc  iea.  with  turctler.cs  and  generally  non  hydro  - 
s  t  a  r  i  c  p  r.  e  r.  c  m  e  r.  a  w  him  otter:  characterize  strong  tie  w  i  r. 
mount  a  incur  terrain  is  procaoly  the  mere  questionable 
aspects  :f  tms  application  cf  the  analysis  technique. 

Oesp i te  tne  limited  success  of  this  application  of 
CCMP1ZX  i".  h. 3. s  *d rf r*  der.cns n'v‘j.u.0ci  to  fc  e  a  useful  6  x  p  e  r  i  rr»  £  n  — 
tal  toil.  Furthermore,  the  economical  utilization  of 
CCMPLF.X  cn  a  small  computer  ( laM  PC -AT!  holds  promise  for 
wider  research  and  operational  applications. 


These  results  indicate  that  the  objective  of  this  study 
has  been  reached,  i.e.,  it  has  been  shown  that  with  appro¬ 
priate  modifications,  a  mass-consistent  wind  model  for  air¬ 
flow  in  complex  terrain  (specifically,  COMPLEX)  can  produce 
realistic  analyses  of  strong  surface  winds  during  severe 
Chinook  windstorms.  Furthermore,  it  appears  that  such  anal¬ 
yses  can  be  accomplished  with  a  limited  observational  net¬ 
work  and  minimal  computer  resources.  The  analysis  scheme 
presented  here  clearly  holds  promise  for  future  operational 
applications.  However,  further  experiments  and  modifica¬ 
tions  are  necessary  to  eliminate  some  uncertainties  in  the 
technique . 

The  shape  of  the  "hydrostatic  lee  wave"  forcing  func¬ 
tion  used  in  the  current  application  did  no:  reproduce  the 
observed  strong  winds  upstream  of  the  lee  wave  trough  axis. 
This  suggests  that  the  flow  surfaces  should  have  been  lower 
in  that  area,  resulting  in  an  overall  configuration  compara¬ 
ble  to  previous  numerical  simulations  of  the  atmospheric 
equivalent  of  a  hydraulic  jump,  rather  than  a  hydrostatic 
lee  wave.  Lacking  detailed  mesoscaie  observations  of  the 
flow  surfaces,  the  use  of  a  two-dimensional,  time  dependent, 
numerical  model  of  the  hydraulic  jump  to  drive  COMPLEX  would 
go  far  in  removing  the  arbitrary  selection  of  the  flow  sur¬ 
face  configuration  (although  microcomputer  resources  may  be 
insufficient  for  such  a  modification) . 


85 


COMPLEX  should  be  further  modified  to  reflect  the  we] 1  • 


mixed  boundary  layer  which  is  characteristic  of  strong  wind 
regimes .  This  change  should  not  allow  the  intersection  of 
the  lowest  flow  surface  and  the  terrain,  thus  eliminating 
the  unrealistic  pattern  of  strong  winds  around  peaks  rather 
than  over  them. 

In  some  of  the  analyses  described  here,  there  was  con¬ 
tamination  of  the  wind  field  at  the  lateral  boundaries  of 
the  domain.  This  problem  can  be  minimized  by  simply  in¬ 
creasing  the  size  of  the  computational  domain  beyond  the 
area  of  interest. 

Finally,  it  is  recommended  that  additional  case  studies 
with  more  observations  should  be  considered  to  determine  (i) 
the  optimum  distribution  of  observation  stations,  (ii)  the 
errors  introduced  by  nor.-hydrostatic  effects,  and  (iii)  the 
relationship  between  the  modeled  winds  and  the  maximum 
gusts  . 
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